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, ^ ■ Abstract 
• 

. The mathematical structure of homogeneous loop quantum cosmology is ana- 

' lyzed, starting with and taking into account the general classification of homoge- 

neous connections not restricted to be Abelian. As a first consequence, it is seen 
qh. that the usual approach of quantizing Abelian models using spaces of functions on 

the Bohr compactification of the real line does not capture all properties of homo- 
bX)' geneous connections. A new, more general quantization is introduced which applies 

to non-Abelian models and, in the Abelian case, can be mapped by an isometric, 
but not unitary, algebra morphism onto common representations making use of the 
' Bohr compactification. Physically, the Bohr compactification of spaces of Abelian 

QQ ■ connections leads to a degeneracy of edge lengths and representations of holonomies. 

. Lifting this degeneracy, the new quantization gives rise to several dynamical prop- 

erties, including lattice refinement seen as a direct consequence of state-dependent 
■ regularizations of the Hamiltonian constraint of loop quantum gravity. The repre- 

, sentation of basic operators — holonomies and fiuxes — can be derived from the 

full theory specialized to lattices. With the new methods of this article, loop quan- 
tum cosmology comes closer to the full theory and is in a better position to produce 
reliable predictions when all quantum effects of the theory are taken into account. 
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1 Introduction 



Loop quantum cosmology [H |2] aims to develop and analyze cosmological models by in- 
corporating crucial guidance from the full theory of loop quantum gravity [HI HI |5]. Even 
though its systems cannot yet be derived completely, they constitute more than a set of 
minisuperspace models. Characteristic effects of quantum geometry have been found in 
this setting, and contact with a potential full framework of quantum gravity has allowed 
one to fix some choices left open in traditional models of quantum cosmology. Nevertheless, 
ambiguities remain in loop quantum cosmology and loop quantum gravity, to be described 
by sufficiently general parameterizations that might be restricted by phenomenological 
analysis. Also for these parameterizations, contact with the full theory is essential: many 
different features collapse on one single parameter when geometry is restricted to exact 
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homogeneity. Disentangling different contributors to one effect is important for estimates 
of typical ranges of parameters. 

Ambiguities notwithstanding, some general features have been found, foremost among 
them quantum hyperbolicity [6t j7j , a mechanism of singularity avoidance based on discrete 
structures of evolution operators. However, effective geometrical pictures of resolution 
mechanisms are difficult to derive since several different quantum effects contribute, in ad- 
dition to higher-curvature corrections also quantum-geometry modifications. (Space-time 
as we know it may not even exist in extreme quantum phases [HI [9].) Specific details of 
physics in deep quantum regimes, for instance the values of upper bounds on the energy 
density of matter, as another possible indicator of singularity resolution, are often unreli- 
able because the setting remains too reduced and too ambiguous. An appropriate viewpoint 
is one akin to effective field theories, where one uses proposals for full (but possibly incom- 
plete) theories to derive generic low-energy properties. With such a view, the framework 
is empirically testable because it can give rise to potentially observable phenomena, even 
if they cannot be predicted with certainty in all their parameterized details. 

Developments in loop quantum cosmology have not always followed a general view, 
especially since the publication of |10], which stimulated a line of research focusing with 
minutest detail on pure minisuperspace models. Strenuous contact with the full theory 
has been replaced by ad-hoc assumptions (for instance related to degrees of freedom and 
scaling behaviors in discrete structures); ambiguities (such as the so-called area gap and its 
postulated dynamical role) have been fixed by hand. Valuable results have been produced, 
chiefly of mathematical interest, showing what discrete features and non-standard quantum 
representations may imply; see e.g. pTl [T2| [13] . For physical statements, however, the 
viewpoint espoused likely contains too many artefacts to be reliable. This article adds 
additional items to the list of known minisuperspace limitations. 

The main aim, however, is to provide a general description of homogeneous reduced 
systems for quantum cosmology, focusing on but not restricted to loop quantum gravity. 
Since quantum cosmological models are beginning to be developed in approaches closely 
related to loop quantum gravity, such as spin foams [HI [HI [16] , it is important to state the 
general setting of quantum cosmological models, and to point out limitations, dangers, and 
promises. The following section begins with a description and classification of symmetric 
models within full (classical or quantum) theories, amended in later sections by specific 
details and discussions within loop quantum cosmology and, briefly, the spin-foam setting. 
The final section will put these models in the general framework of effective theory. Along 
the road, we will be led to several mathematical features overlooked so far. Most impor- 
tantly, Hilbert-space representations based on functions on the Bohr compactification of 
the real line do not properly capture all aspects of homogeneous connections; rather, they 
give rise to a degeneracy of two important parameters corresponding, in the full theory, 
to the edge length and representation label of holonomies. The origin of the degeneracy 
is identified here as a mathematical coincidence realized in Abelian models only. To solve 
these problems, we will present a new non-Abelian construction, work out a detailed re- 
lation to the full theory, and arrive at a new viewpoint on dynamical difference equations 
with a natural implementation of lattice refinement. 
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2 Reduction 



A classical reduced model, realizing a given symmetry, uses an embedding Ai ^ S from 
the set Ai of symmetric geometries into full superspace S. A minisuper space geometry 
as an element of Ai is specified by finitely many parameters aj (which, to be specific, 
one may think of as the three scale factors of a diagonal Bianchi model), mapped to a full 
metric gab{(^i) by the classification of invariant metric tensors. (A classification of invariant 
connections and triads is mathematically more well-developed [T71 HB]; see [ISl EDj for 
applications in the present context.) Inserting gabiaj) in equations of motion, the action or 
constraints of the full theory then produces corresponding equations for the finitely many 

We distinguish between a minisuperspace model (a definition of minisuperspace A4 
with a dynamical flow on it) and the stronger notion of a reduced model (a minisuperspace 
model including also an embedding Ai ^ S, making contact with the full theory). There 
is not much of a difference classically, but there is a big one at the quantum level. In 
quantum cosmology, most models remain in the minisuperspace setting, defining some 
dynamical flow on a system with finitely many degrees of freedom such that the dynamics 
of general relativity follows in a semi classical limit, perhaps with inspiration by but no 
derivation from some full theory that quantizes S. The key ingredient of reductions — 
making contact with a full theory of quantum gravity, if only in a weak sense — most often 
is missing. This article deals with the problem of reduction at the quantum level, going 
beyond pure minisuperspace models. 

2.1 Reduction, selection, projection 

In addition to minisuperspace versus reduced models, it will be useful to distinguish be- 
tween three classical procedures of deriving symmetric systems beyond a mere minisuper- 
space prescription: 

Reduction: The definition of a reduced model, as already stated, contains two parts: 
embedding in 5 and deriving a dynamical flow by inserting minisuperspace metrics 
into the full equations. After a reduced model has been defined, one can proceed to 
solve its equations and evaluate solutions for potential predictions!! 

Selection: Instead of starting with a set of minisuperspace geometries and deriving a 
reduced flow from the full theory, one could first solve the full equations and then 
select symmetric ones that admit a set of Killing vector fields with an algebra of 



^The procedure of inserting symmetric tensors may not commute with variations used to compute 
equations of motion [2TJ[22]: symmetric actions or constraints do not always produce the correct symmetric 
equations of motion. In our following discussions we make use of constraints, and therefore must assume 
that variation and symmetrization commute. In the Bianchi classification of homogeneous models, for 
instance, we are restricted to class A [53]. 

^ "Embedding" might be a better name for this prescription, but "reduction" is much more standard. 
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the desired symmetry type. In general relativity with its complicated and largely 
unknown solution space, this procedure is rather impractical. 

Projection: Or, again starting with solutions to the full equations, one may derive a 
symmetric geometry for every full solution by some kind of averaging process. In 
addition to the problem of a selection procedure, one would have to face the daunting 
averaging problem well-known from cosmology fM\ . 

While reduction is a standard classical procedure, selection or projection cannot be per- 
formed by current means. In quantum cosmologyjfl the best one may attempt is therefore 
a quantum version of reduction, going beyond pure minisuperspace quantizations but not 
directly accessing the full solution space. This is the topic of the present article. 

2.2 Quantum cosmology 

Several difficulties arise when one tries to extend classical reductions to quantum theory. 
Unlike classically, symmetric solutions can no longer be exact: Inhomogeneous degrees of 
freedom are set to zero, both for configuration and momentum variables. The spatial metric 
at any time must be invariant under the given symmetry, and so do its rate of change or 
extrinsic curvature for the metric to be able to remain invariant. Setting non-symmetric 
modes of both canonical fields to zero violates the uncertainty principle, and symmetric 
quantum solutions cannot be exact solutions of the full theory. Reduced quantum models 
can at best be approximations, but making sense of the approximation scheme in a clear 
way remains one of the outstanding challenges to be faced. 

Mimicking classical constructions, the kinematical structure of quantum gravity can be 
reduced by making use of a mapping a: "Hhom — ^ ^^fuii from the kinematical Hilbert space 
"Hhom of a homogeneous minisuperspace model, quantizing the degrees of freedom a/, to 
the space of distributional states in the full theory |19j. This mapping is analogous to 
the classical A4 ^ S, but the distributional nature of the target space spells additional 
complications. Moreover, with uncertainty relations violated, symmetric quantum evolu- 
tion is not exact in the full theory. Starting with a homogeneous full (but distributional) 
state V'hom G o"(^hom) C I^fuu, the distributional extension of (the dual action of) the full 
constraints Cfuii or their gauge fiows exp(i(5Cfuii) does not leave the state in the image of 
homogeneous states. In loop quantum cosmology, methods exist to define and analyze 
maps (70 

but the derivation of a dynamical fiow from the full theory remains difficult even though 
candidates do exist. Loop quantum cosmology therefore realizes an incomplete quantum 
reduction. Wheeler-DeWitt quantum cosmology, on the other hand, is not a reduction but 

■^A great danger in quantum cosmology is that its procedures amount to neither reduction nor selection 
nor projection, but rather to production — not a combination of projection and reduction as the word 
might suggest, but just the presentation of an artificial model of unknown pedigree. 

"'More precisely, as detailed below, instead of Pfuii a distributional space based on lattices, and therefore 
fully inhomogeneous but not the most general states, are used. Otherwise, if edges not adapted to the 
symmetry appear, there are obstructions to embeddings of states |26| . 
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a pure minisuperspace quantization since no analog of a exists. (There are also models of 
loop quantum cosmology which do not go beyond pure minisuperspace models, disregarding 
proper considerations of a.) 

In order to restrict or truncate full quantum evolution (or gauge flows) to homogeneous 
states, one must specify a projection of CfuuV'hom or exp(i5C*fuii)'?/'hom back to the image of 
Hhom in ^^fuih for all states V'hom and all full constraints — some part of the averaging prob- 
lem plays a role even for reduction when quantum effects are involved, another indication 
of more-severe problems. No such projection has been provided so far, and therefore the 
dynamics of reduced models, let alone minisuperspace models, remains incomplete. 

The problem is challenging not just owing to quantum issues, such as the distribu- 
tional nature of symmetric states. Even classically, the question of how to project a non- 
symmetric metric to a homogeneous one is complicated, and unresolved in its generality; 
it constitutes the averaging problem of cosmology. Since a complete derivation of reduced 
quantum models from the full theory of some form would, in its semi classical limit, include 
a solution to the averaging problem, one cannot expect progress on the dynamical front of 
quantum reduction unless the classical averaging problem is better understood. 

The averaging problem remaining open, the only way at present to go beyond minisu- 
perspace models is to use properties of a homogeneous background for an implementation 
of inhomogeneity, perhaps by perturbation theory. In classical cosmology, one commonly 
makes use of this perspective when inhomogeneous fields are expanded by Fourier transfor- 
mation with respect to the modes on a homogeneous background. Classically, the approx- 
imation is well-understood. In quantum cosmology, the procedure suffers from the same 
problems encountered for homogeneous models, and adds new ones related to the quanti- 
zation of inhomogeneous modes. Also the question whether results may depend sensitively 
on the background (and often gauge) chosen before quantization remains thorny, related to 
the complicated anomaly issue of quantum gravity. (Some anomaly-free realizations with 
partial quantum effects exist [271 12H1 ISll [30l [31], [32], [33] . The anomaly problem is to be 
faced in canonical and covariant approaches alike: in canonical approaches it appears in 
the constraint algebra; in covariant ones, in the path-integral measure or in the correct 
choice of face amplitudes of spin foams [M]-) 

Facing these difficulties, it is an effective viewpoint which allows progress, making use 
of sufficiently general parameterizations of quantum effects, but disregarding fine details. 
The viewpoint is half-way between minisuperspace models and a complete dynamical em- 
bedding in the full theory: One avoids the averaging problem by using inhomogeneous 
model states adapted to the symmetry, much like classical Fourier modes on a given back- 
ground. In practice, it is often lattice states with links along the symmetry generators 
that allow one to include a discrete version of inhomogeneity at the quantum level [35] . A 
background structure is then built into the framework, but it becomes possible to deal with 
inhomogeneity, going beyond pure minisuperspace models and escaping their limitations 
and artefacts. A background or some gauge fixing has entered, possibly giving rise to new 
problems. But at this stage, effects sufficiently general and parameterized can give reliable 
access to the physics involved. 

In loop quantum cosmology, this procedure has been developed to the degree that 
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cosmological phenomenology can be done. The kinematical structure — the basic algebra 
of holonomy and flux operators — can be derived from the full one. Evolution and the 
dynamics, facing the classical averaging problem and the anomaly problem of quantum 
gravity, remain much less understood, but here parameterizations have been developed 
that capture interesting effects. Especially the interplay of various quantum corrections, 
signiflcant in different regimes of curvature, puts restrictions on possible phenomena. 

We will now go back to the basics of these constructions to clarify and generalize 
several mathematical objects involved. We will point out one major problem due to oft- 
used Abelianizations of cosmological models, overlooked so far. Its solution has several 
ramiflcations even at the level of formulating the dynamics. 



3 Loop quantum cosmology 

Loop quantum gravity represents the basic geometric flelds, connections A^^ and densi- 
tized triads , by integrated versions: holonomies of the connection and fluxes of the 
densitized triad. As a flrst step toward symmetry reduction, mathematical theorems are 
available to classify and construct different types of connections invariant under the ac- 
tion of some symmetry group [17J. When specialized to homogeneity [36], or a transitive 
group action on space, a set of models equivalent to the usual Bianchi classiflcation re- 
sults: For every Bianchi type, there is a set of three left-invariant 1-forms cuf , 7 = 1,2, 3, 
such that invariant connections are = c}a;f with spatial constants (but time- dependent) 
c}. Invariant densitized triads take the dual form, E°- = p[X"|deta;^| with invariant 
vector flelds Xf dual to cuf. This choice of densitized-triad components ensures that 
{87T'~fG)^^ Al^E^d^x = {8'n-'~fG)^^Voc\pl, integrating the symplectic term of an action 
over some region V of coordinate volume Vq = /y d'^x. Up to constant factors, c} and pf 
are therefore canonically conjugate: 

{a;.p/} = ^*]*/. (1) 

Vo 

A Bianchi metric of the given type is qab = \ det{p^)\p}pjU^LO}l with inverse matrices 
of p[. The metric is invariant under rotations pf i— > Rlpj with R G S0(3), generated as 
gauge transformations by the Gauss constraint eij^c'jpl^. The diffeomorphism constraint is 
not relevant for homogeneous models, and the Hamiltonian constraint is [37] 

8nG^\det{pl) \ V 

+2(1 + 7-^)(£i - T^M - T'})P^P^) 

with the structure constants Cfj of the Bianchi group and the spin connection F}, depend- 
ing on Cfj and p^. In diagonal models, where c} = C(/)(5} and = p^^^S^ (no summation 
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over J), ([2]) reads [38] 



H = ^(((c2r3 + c3r2-r2r3)(i + 7"')-riV-7"'c2C3)/bVM 

+ {{c,r, + C3T^ - r,r,){l + 7-^) - nh^ - 7~'ciC3) V¥WM 

+ {{c,T2 + - rir2)(i + 7'') - n'cs - r^c,c2) VIpVM) (3) 

with spin-connection components 

and coefficients of the Bianchi classification. 

For Bianchi I with Cfj = 0, the constraint reduces to 

_ 1 c>iCjPjPi - c%plpi 

Sn^'G ^\det{pi)\ ■ ^ ^ 

In a spatially fiat isotropic model with c] = c6} and pf = p6- with c = 'ya and |p| = 

— see [361 [I] fo^' the origin of the factor 1/4 — this expression reduces correctly to the 

gravitational contribution H = —3{A7T'y'^G)^^c? ^/\p\ of the Friedmann equation. 



3.1 Abelian artefacts 

One of the first steps of loop quantization consists in replacing connection components c} 
with holonomies or exponentials exp(c}rj) G SU(2). The vast majority of investigations in 
homogeneous loop quantum cosmology, however, deals with Abelian models in which the 
original SU(2) is replaced by U(l), thanks either to additional isotropy symmetries [39j or 
a diagonalization assumption [38]. With isotropy or diagonalization, a classical invariant 
connection automatically becomes Abelian, and the use of U(l) is not ad hoc but required. 
However, in quantum cosmology, Abelian structures turn out to allow specific choices of 
Hilbert-space representations not possible in non-Abelian ones. Such quantizations, based 
essentially on spaces of functions on the Bohr compactification of the real line, are therefore 
in danger of introducing additional artefacts, structural properties that cannot be met in 
non-Abelian models, let alone the full theory. 

For instance, an isotropic connection has the form A]^ = c6l^ with just one phase-space 
component c [39]. Mimicking matrix elements of holonomies along straight lines such as 
the edges of the integration cube V (or some other fixed set of edges), one first represents c 
by U(l)-holonomies h = exp{iVQ^^c), or = exp(mc) with an integer U(l)-representation 
label n and c := Vq^^c. By spanning a function space with superpositions of /i" for all 
integer n, all continuous functions on the group U(l) are realized. 

From U(l)-holonomies one can reconstruct the connection component c only modulo 
27r. One gains full control over the connnection if one considers holonomies along all pieces 
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1 /3 

of the edges of the integration cube of lengths < K) 5 such that holonomies = 
exp(i/ic) with /i G M result, where /i may be considered as a product \n of the fractional 

1 /3 

edge length A = ^q/Vq with the representation label n. Allowing for superpositions of 
all h'^ as an orthonormal basis, the Hilbert space of all integrable functions on the Bohr 
compactification Meohr of the real line is obtained |10], rather than a function space on 
some periodification of M. 

In this procedure, which has become standard, one implicitly makes use of an identity 
realized for representations of Meohr but lacking a non-Abelian analog. In the Abelian 
case, we start with the U(l)-holonomy exp(2Ac) and evaluate it in the p„-representation: 
Pn{exp{iXc)) = exp{iXnc). It so happens that this is the same function of c as obtained 
from pAri(exp(ic)), now using a representation of Meohr- Holonomies in the n-representation 
have led us to the first expression, which is then identified with the latter. Since they agree 
as functions, one may base Hilbert-space constructions on functions on Meohr- However, 
this step is not available for non-Abelian models, in which case there is no relationship 
between pj{exp{\A)) and pxj{exp{A)) for A in the Lie algebra of the group, usually SU(2). 
The second expression is not even defined unless A is an integer (or a half-integer if j is 
integer), but even then, the two matrices are unrelated. 

If functions on Meohr are used, one must proceed with care to avoid artefacts in Abelian 
models, a problem which has not been realized in existing constructions. Mathematically, 
one would confuse p„(exp(iAc)) with pxn{exp{ic)), which are identical as functions of c but 
have different meanings and are elements of different function spaces. Physically, merging 
A and n to one number p = Xn, as done when the Bohr compactification is used, eliminates 
important information because the edge length A and representation label (or geometrical 
excitation level) n are then indistinguishable. In operators, however, A and n should play 
rather different roles according to what is known from the full theory. In this section, we 
present a new quantization of homogeneous models in which A and n are kept separate, 
lifting their degeneracy, and which extends to non-Abelian models. 

3.2 Homogeneous holonomies 

A local homogeneous connection Ahom is a 1-form on the translational symmetry group S 
underlying some Bianchi model, taking values in the Lie algebra CG of the structure group 
G, G = SU(2) for gravity in Ashtekar-Barbero variables. If the symmetry group acts 
freely, without any isotropy subgroups, there is a one-to-one correspondence |17j between 
homogeneous connections according to S and linear maps 0: CS — )■ CG (or elements of 
CS* X CG), not required to be Lie algebra homomorphisms. Given (p, the corresponding 
homogeneous connection is the pull-back Ahom = (p*^MC under (p of the Maurer-Cartan 
form, which latter can be written as cjmc = oo^Tj in terms of left-invariant 1-forms = 
cjfdx" on 5* and its generators Tj. The homogeneous connection components introduced 
before are the coefficients in 4>(Tj) = c]Ti with generators Tj of CG (here, Tj = —\icFj in 
terms of Pauli matrices) . 

A minisuperspace model quantizes the components c}, or rather the linear maps cf). Both 
ingredients are in one-to-one correspondence, but the additional structure shown by the 
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linear maps is useful to decide how different quantum numbers, such as A and n in Sec. 13.11 
should be related to properties of S (space) and G (internal space). We will therefore 
derive a quantization based on the mathematical structure of (f). To extract independent 
degrees of freedom, we fix a set of generators T/ of CS and understand a homogeneous 
G-connection for a given syrnmetry group 5* as a set of maps 0/: (T/) — )■ CG with the 
scaling condition (j)i{rX) = r(f)j{X) for all r G M. 

Following the methods of loop quantum gravity, we quantize connection components in 
terms of holonomies. According to the structure of homogeneous connections, we introduce 
the notion of homogeneous holonomies by exponentiation — maps hfp: CS — )■ G, h^{X) = 
exp(0(X)) with the scaling condition h^{rX) = hr^{X). The maps 0/ = used 
here differ from 0/ by factors of L/, side lengths of the integration region V of volume 
Vq = L1L2L3, assumed cubic (spanned by the three generators Tj = Xfd/dx"" of the S- 
action). If the sides of V are aligned with the three symmetry generators, one may think 
of h^{Ti) as the holonomy along the corresponding side. This relationship will be made 
more precise below. 

Elements X G CS to which is applied carry information about the edge used to 
compute holonomies h^{X). Referring to the Killing metric on CS, we decompose X = 
Xv into its norm A = |X| and the unit vector v = X/\X\, corresponding to the 
coordinate length of the edge and its direction. With the scaling condition, we then have 
h^{X) = exp(A0(f)). We can compute all information about from derivatives 0(Tf) = 
d/i<^(AT7-)/dA|A=o- We are indeed representing the space of all homogeneous connections, 
not some periodic identification. 

The dependence of homogeneous holonomies on A = |X| will play an important role in 
our constructions. If we consider an edge e/ of coordinate length ij along the generator 
Xj, the holonomy, in general a path-ordered exponential he = V exp{J^dsAl^Tie°') of the 
connection integrated over a spatial curve e(s), is h^j = exp(£/c}rj) = exp{£jLj^ (f)(Tj)) = 
h(f,{XjTi) with A/ = ii/Lj. If all edges are contained in the integration region, we always 
have A/ < 1. More generally, we can allow all real values, but for SU(2), given periodicity 
of the exponential function, may restrict to < A/ < 47r without loss of generality. 

For a given connection, the three choices for /, or three directions of space, give rise 
to three independent SU(2)-elements h^{Tj). For fixed A/, one can therefore describe the 
space of homogeneous connections in terms of SU(2)^jfl but the connection used can be 
reconstructed completely from the holonomies only if different choices for A/, or curves 
of different lengths, are considered. If the curves and their lengths are fixed, as in the 
original constructions of [3Sl ESI ES], only a certain periodic identification of the space of 
connections is realized. 

Dropping the reference to particular edges, minisuperspace models refer to the following 

^Thanks to homogeneity, each holonomy transforms as h^j i~> gh^^g^^ under an internal gauge trans- 
formation, with the same g € SU(2) for all three edges and on all their endpoints. These transformations 
are identical to those obtained in the full theory for three closed loops intersecting in one 6-valent vertex 
|36) . One may picture homogeneous spin- network states as such vertices, but with homogeneity, the vertex 
corresponds to all of space — homogeneous states are distributional and not given by single spin networks; 
see SecEXl 
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structure: Our homogeneous holonomies h^{X) are elements of a space of functions gf. M x 
CG — i- G, {Li, (pi) h-). gi{Li, (pi) required to satisfy giijLi, r~^(j)i) = gi{Li, (pi) for all r G R 
(the scaling condition). The choice r = 1/L/ shows that any such function can be written 
as gi{Li,(pi) = gi{Li(pi) with a function gi of just one variable A G CG. If L/ and r 
are fixed, gi is simply the group exponential; setting r free allows for different scalings or 
different sizes Li of the integration region within one model. 



3.3 Representation 

For homogeneous models of loop quantum cosmology, we turn the function space based on 
holonomies into a Hilbert space with an action of holonomies and fluxes as basic operators, 
such that their commutator corresponds to the classical Poisson bracket ([1]). One immedi- 
ate problem caused by the non-Abelian nature of general connections in combination with 
homogeneity regards the way of exponentiating connection components to holonomies and 
obtaining a closed basic algebra for {exp(A7c}ri), p^}. Once the path-ordering of inhomoge- 
neous holonomies is no longer available, derivatives of exp(A/c}rj) by Cj will produce extra 
factors of Tj between products of c}rj in a power-series expansion of the matrix exponential: 



J n=0 k=l 



4r,)^ ^Tj{c\n) 



but they do not automatically factor into products of exponentials with Tj to mimic the full 
holonomy-fiux algebra. (While the cotangent bundle T*G defines a natural phase space 
with group- valued configuration variables, it does not necessarily model the correct relation 
to inhomogeneous holonomies.) 

For a closed basic algebra to result, the factors in derivations of basic holonomy-like 
functions of c} may have to be re-ordered, but within a pure minisuperspace model, there 
is no guideline, no trace of the path-ordering left by which one could construct a natural 
ordering. By looking more closely at the relation between basic operators in models and 
the full theory, we will be led to one distinguished choice. 



3.3.1 Function space 

We first construct the space of all functions on homogeneous connections gi and turn it 
into a C*-algebra. According to the definition of gi, we look at functions ip of L/ and 0/ 
such that the dependence happens only via gi{Li,(pi). The scaling condition for gi then 
translates into an analogous condition for ip. 

If we fix Li, considering gi{Li,(pi) simply as an element of G, and refer to the Peter- 
Weyl theorem, the general dependence on 0/ can be realized by superpositions of functions 
{m\pj{gi{Li, (j)i))\n) with all irreducible representations pj of G and elements |m) and \n) 
of an orthonormal basis of the representation space of pj. Setting Li free, with gi{-, (pi) as 
a 1-parameter family of G-elements, a larger class of functions is possible. However, the 
scaling condition can be realized only if our functions are superpositions of Pxj{gi)n '■ = 
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{m\pj{gi{\Li,(j)i))\n) for A G Q. (The restriction to rational as opposed to real A will be 
motivated momentarily. The labels A, j, m and n may depend on J, but we will often 
suppress the dependence for notational simplicity.) 

We multiply two functions Pxijiigi)^^ and P\2,j2i.9j)7i2 follows: If / 7^ J, we simply 
take the product function depending on gj and gj as independent variables, thereby gen- 
erating a tensor-product space. If / = J, we write Ai = p/r and X2 = q/r with the least 
common denominator r of Ai and A2 and define 

Px„n{gi)Z' ■ P>^2,hi9i)n2 ■■= Yl PyrjA9i)Z'Pyr,ni9i)hl---pi/rjA9itr" 

hi,...,hp^l,ki,...kq^i 

^Pl/rj2i9l)Z'Pyr,ni9l)ll ■ ■ ■ Pl/r,j2i9l)nV ' (6) 

One may decompose the products of matrix elements on the right-hand side into superpo- 
sitions of irreducible contributions to the tensor product pf^ ® pf^, akin to a spin-network 
decomposition in the full theory. The product is then again a superposition of Pxjigi)^- 
Multiplication as defined is commutative and associative because these properties are re- 
spected by the least common denominator. There is a unit element given by A = (in which 
case the value of j does not matter). For vanishing j, having the trivial representation of 
SU(2), one may expect a trivial action, too. However, according to IQ, multiplication with 
Pxfiidi) for A 7^ may still give rise to decompositions of factors, providing a different form 
of the function product even though the values taken by the original function and its prod- 
uct with pxfi{gi) do not differ. The functions pxfligi) (matrix indices are not required in 
the trivial representation) play the role of refinement operators, decomposing a holonomy 
into pieces whose length is determined by A in relation to the corresponding parameter 
of the state acted on. Since pQ{h^{XTj)) = 1 classically, these refinement operators have 
no classical analog, as one may expect for a classical theory knowing nothing about the 
underlying discreteness. We will soon make good use of pxfligi)- 

The multiplication rule observes homogeneity: One can interpret the law as a decom- 
position of two initial holonomies of different lengths as products of equal-length pieces. 
Without homogeneity, these pieces at different places would be independent, but homo- 
geneity makes them identical. We therefore take the tensor product of all small pieces, split 
as illustrated in Fig. [H and sum over indices according to the product form. At this stage, 
it becomes clear why A should be rational: For incommensurate Ai and A2, the product 
PAiji(s/)™i^ ■ Px2,j2i9i)ri2 would have to be split into infinitely many infinitesimally small 
piecesO 

We define a star operation by (pxjigi)^)* '■= Pxjigi)^ (related to matrix elements of 
the dual representation p* of pj for unitary groups). The space of functions turns into an 
Abelian C*-algebra with the supremum norm, assuming G to be compact. 

A Hilbert-space structure is obtained by combining the product rule with the Haar 

^Projective-limit constructions, in which states with a given denominator r would play the analog of 
fixed-state cylindrical states in the full theory, may be used to allow for incommensurate parameters A, 
but we will not need this in the present article. 
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Figure 1: Evaluated in homogeneous connections, different pieces of equal length in one 
holonomy amount to the same function (left). To avoid overcounting degrees of freedom, 
multiplying a holonomy with a shorter one requires a decomposition into pieces of maximal 
common length (right), giving rise to the rule 



measure on G. We define the inner product 

where r is again the least common denominator of Ai and A2. On the right-hand side of 
this equation, the value of r no longer matters because we integrate over all group elements 
gKiLx/r, 0x), the sole arguments of ipx^jAgi)^^)* ■ Px2,j2i9.j)Z^ ■ 

3.3.2 Fluxes 

Components pf = LjLkPI {^ijk = 1) of the densitized triad, canonically conjugate to 
c} = L/c} via {c},p^} = SiTjGdjd^j, are quantized to operators with action 

plpxA9j)n ■= -8vr274 >^S'jPxAr^9j)n (8) 

where we define the short form p\,j{'rigj)n ■= Zlfc Pi(^i)™PA,j(5'j)n- 

Non-Abelian flux operators as defined are not symmetric because the product ([6]) in ([7]) 
in general includes a decomposition. (One may think of the measure factor as including 
refinement operators Pi/r,o{9K)-) While a flux operator acting on either entry in the inner 
product inserts a Xj to the left of gj according to ([H]), the integration required to evaluate the 
inner product splits holonomies according to one divided by the least common denominator 
of Ai and A2. Integration by parts, performed after multiplying according to and 
decomposing, would then insert in each decomposed contribution, not just to the left 
of the whole gx- Also here, projective-limit constructions as already alluded to can be 
used to define self-adjoint non-Abelian flux operators for homogeneous models, but we 
refrain from working out the details in this article in which we are interested primarily in 
the relation between models and the full theory. Non-symmetric fluxes then arise because 
reduction entails averaging and decomposition, and what appears as a simple quantization 
of the densitized triad in a pure minisuperspace model turns out to be a more complicated 
operator when inhomogeneous degrees of freedom are taken into account. We will see this 
in more detail in Sec. 13.4.31 Even non-symmetric flux operators (quantizing the densitized 
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triad combined with the action of holonomy- dependent decomposition, reahzed for instance 
by Pxfiigi)) model the behavior of the full theory, in which fluxes are self-adjoint. 

Also the commutator [pA,j(fl'j)™,pf] of basic operators in the non-Abelian holonomy-flux 
algebra requires care due to decomposition. Up to ordering, it equals 8'7ii'yipXSjPxj{Tigj)'^ 
and corresponds to ih times the Poisson bracket of pj{exp{XcjTi))^ , the classical analog 
of px,j{gj)ny s-iid pj. For an example of ordering issues, as indicated in the beginning of 
Sec. 13.31 look at the commutator [pi/2,j(fl'/)jv ^pf] acting on the state ip{gi) = Pi,j{gi)n- 
Acting with pj first produces a single insertion of Tj to the left of gj. Acting with pf after 
the action of pi/2,j{,gi) produces two insertions because we first decompose pi/2j(fi'/)™ ■ 
Pi,jigi)n = T^k Pi/2,jigi)n'' Pi/2,jigi)TPi/2,jigi)n according to dHD and use the Leibniz rule. 
The result, 



A7ri-filpi/2,jiTigi)^ilj{gi) + ATri-filpi/2jigi)N ^ pi/2,j{gj)k Pin^Mgifn - Piji'^igi] 



m 



is as expected up to the second term, a contribution that can be made to vanish by 
reordering. 

Similar calculations for arbitrary Ai and A2 in [pAiji(fi'/)jv ?Pf] acting on a state ipi^gi) = 
P\2,j2i.gi)n (^-contributions with J 7^ / do not matter) result in 

[PxMn.pI] = 87,il(l\px,j{ngi)N + ^thi^Ib'.pxMn ■ (9) 

The specific ordering introduces an extra contribution from the reordering operator Rj 
defined by 



1 ^ 

Ri{pxi,hi9i)Z pxN,jMigi)Z) = -^p^iAgi 



n=l 



Yl pyr,Mk: ■ ■ ■ pi/r,jAngi)ll^, ■ ■ ■ pxir,,Sgi)lY-' 

, p=l A;i,...,fcA„r-l 

■ P^NdNigi) 



N 



njv 



n=l 



with the least common denominator r of all A„ involved. Up to ordering, -Rf (paiji(5'/)™^ ■ 
■■■ ■ P>^N,jNi.gi)n^N^ vanishes. For more compact notation, one may express the refine- 
ment included in the action of R{ in terms of pxfl{gi) with suitable A. Therefore, i?f 
is not independent of the operators already introduced. Specifically, we can write Rj = 
— {S7ii'yip)~^pl{pi/rfl{gi) — 1) where 1/r is determined as the least common denominator 
of all A-parameters of holonomy factors to the right of Pi/r,o{gi)- The basic commutator 
then reads 

[pxAgi)N,pi] = -ipipxAgi)N)- Piipi/rM - i)PxAgi)N = {p,ih^{mrj,Pi}- m 
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Since p\fl{gi) = 1 classically, the commutator can indeed play the role of a quantization 

of the Poisson bracket, {pj{h^{XTj))^ ,pj} = SnijipXpxjirigf)^ . The basic algebra is 
more complicated than one may have expected, but we do obtain a closed algebra of basic 
operators with the correct classical limit. 

Each operator pj can be viewed as a component of the flux operator for a surface given 
by the full side of a cubic V, transversal to the direction Xf. These minisuperspace fluxes 
refer just to the artificial integration region V and its coordinate volume Vq = L1L2L3, not 
to actual edge lengths or areas of dual surfaces, for instance in a lattice. One can easily 
rescale the linear flux by using \j\kp{ = ^jd-Kpl for ^ijk = 1, now corresponding to the 
flux through a surface of area related to edge lengths £j = \jLj and Ik = ^rLr- Without 
a lattice construction to show how dual surfaces are related to links, however, a strict 
minisuperspace quantization does not provide a satisfactory, V-independent definition of 
fluxes. We will complete the construction in Sec. 13.41 when considering the relation to the 
full theory. 

A convenient set of states in the homogeneous Hilbert space is given by a form of 
spin-network functions, depending on the connection via finitely many holonomies. 



(fe) 



(11) 

with coefficients il).(k) ,(k) .(fe) .(k) .(/=) .(k) (k) (k) (k) (k) (k) (k) for k in some finite index set. 

Aj jAg ,J-^ ,J2 J3 ,"^1 ,"^2 '™3 '^1 '^2 '"3 

Acting on the contribution Pxjjj{gi), Jl = (8vr7£|)^^A7^pf satisfies the su(2)-algebra by 
the definition in (|H]). The fiux spectrum is therefore given by all numbers Sn'yipXm with 
half- integer m, and the area spectrum (or the spectrum of \/p(p^) by 8^T'yipX^/j{J^^l) ■ 
Although these eigenvalues are real, one can see the non-symmetric nature of non-Abelian 
fiux operators: Eigenstates with different eigenvalues (specifically, states with the same j 
but different A) are not necessarily orthogonal, again owing to the decomposition in ([7]). 
With all rational A allowed, these spectra form continuous sets, but all eigenstates are 
normalizable. The spectra are pure point. 



3.3.3 Properties 

The construction has several convenient features: 

• Because is required to be a linear map, we need consider holonomies only along lin- 
ear combinations of the generators Tj, identified with tangent vectors along "straight 
lines." No path-ordering is required to compute these holonomies, and they have 
simple forms. (For examples of more complicated holonomies, see [26].) In particu- 
lar, their matrix elements span a non-Abelian version of the space of almost-periodic 
functions: superpositions of periodic functions with different periodicities given by 
A. 
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• The space of homogeneous holonomies is compact if G is compact, as the spectrum 
of a unital Abehan C*-algebra. Homogeneous loop quantum cosmology makes use of 
function spaces on compactifications of the classical spaces of homogeneous connec- 
tions, which contain all classical connections as a dense subset. 

• If G is Abelian, the homogeneous Hilbert space is (non-unitarily) related to a product 
of the Hilbert space "Heohr of square- integrable functions on the Bohr compactification 
of the real line, with dim(S') x dim(G) factors. More details will be given below. 

• If 5 does not act freely, its isotropy subgroup requires additional identifications of the 
components of Linear maps must then satisfy </)oad/ = adi?(/)0 for any / in the 
isotropy subgroup and a corresponding element F{f) G G; see [HI US ED] for details. 
Accordingly, we restrict homogeneous holonomies by h^^&dfX) = /iad^(j)O0(-^) in 
addition to the scaling condition. These functions, and therefore gj, take values in a 
subgroup of G, the centralizer of the subset of all F{f) in G, which is Abelian if the 
isotropy subgroup is sufficiently largeJil 



3.3.4 Diagonalization 

Detailed ordering prescriptions make some formulas of representations of non- Abelian mod- 
els look rather complicated. It is more straightforward to describe the space of connec- 
tions by holonomies when diagonal Bianchi models are used, implying Abelianization. In 
diagonal models, the relevant parameters appear by writing a homogeneous connection 
as Al^ = C(/)a;fA} with left-invariant 1-forms cjf of the given symmetry type and S0(3)- 
matrices A} which can easily be fixed to equal 6}. (Writing c] = C(/)A}, not summing over /, 
makes use of the polar decomposition of matrices |38].) The conjugate field, the densitized 
triad, has a similar decomposition, = p^^^XfAj \ det uj^\. All ingredients except cj and 
are determined by the symmetry type or gauge choices, and cj and are the canonical de- 
grees of freedom. As before, the symplectic term (87r7G)~^ (PxAl^Ef = Vo{8'7i'yG)~^cip^ , 
integrated over some bounded region V of volume Vq, provides Poisson brackets 

lc,.f}-'-^S{. (12) 

Holonomies hej = exp(^7C(7-)A}rj) of diagonal connections still take values in SU(2). 
However, the relations between h^j for different / are not arbitrary because their generators 
Yj := ijC(^i)A}Ti satisfy the condition tT{YiYj) = for / 7^ J. (Different rows or columns 
of the S0(3)-matrix A} are orthogonal with respect to the SU(2)-Killing metric rjij = 
— 2tr(rjrj) = 6ij.) Although the h^j do not commute with one another, any pair of them 
obeys gh = hg + h~^g + hg^^ — tT{hg) [SB]: A product of diagonal holonomies with h 
appearing to the right of g can be reordered so that h appears on the left in all terms. The 
gauge structure of diagonal models is essentially Abelian. 



^Formally, classical Abelianization may also be implemented via second-class constraints [HI 142) . 
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To make Abelianization manifest, one usually works with matrix elements hj = exp{iiiCi/2) G 
U(l), completing the reduction of the theory to an Abelian one. For diagonal mod- 
els, all phase-space information is indeed captured by these matrix elements because 
exp{iiC(i)A]Ti) = cos(£/C//2) +2A}rj sin(£7C(/)/2). Fluxes computed for surfaces normal to 
X'^ are = ijijp^ {^i.jk = !)• Holonomies as multiplication operators on states in the 
connection representation and fluxes as derivatives simplify significantly by Abelianization. 
For instance, the volume operator, notorious for its complicated spectrum in the full theory 
[13| HI] and also on the 6-valent vertices of non-diagonal homogeneous models [15] , is a sim- 
ple product V = \j \F^F'^F^\ of derivative operators F^ = —Sni'yEpXjXjd/dcK {^ijk = 1) 
on U(l). Abelianization also implies that a triad representation becomes available. 



3.3.5 Abelian homogeneous connections 

In diagonalized or isotropic models we encounter Abelian homogeneous connections with 
G = U(l). In this case, the structures introduced for general homogeneous connections 
simplify: Our function space consists of superpositions of functions px,n{,g) of a single 
variable g{L,(j)) = exp(iL0) per independent direction, with A G Q and n G N, using 
U(l)-representations Pn{g) = 9^- Multiplication ([6]) now reads 

PAi,ni {g) ■ P\2,n2 (g) = Pl/r,ni {gTPl/r,n2 {gY = Pl/r,pm+qn2 i (13) 

again with the least common denominator r of Ai = p/r and A2 = q/r. The star relation 
is p\^n{g)* = Px-n{g), and the inner product (JTj) evaluates to 



ipx^,mig),Px2,n2ig)) = TT d(0/r)pAi,-ni(exp(i0)) ■ pA2,n2(exp(#)) (14) 

2vr Jo 

— / axp-pni+qn2\^^V\'^^ ) ) — 0pni,qn2 — S i \ „ _ \ „ 

ZTT Jo I ^ MTT-l — A2^2 

after substitution. Finally, the derivative operator (jH]) is 

PPxAa) = Sn-fEl Xnpx,nig) (15) 

with eigenvalues SnjipXri. 

These equations bear some semblance with representations on function spaces on the 
Bohr compactification of the real line, but they are not identical. 



3.3.6 Relation to the Bohr compactification of the real line 

As recalled in Sec. 13. H traditional loop-based minisuperspace quantizations, for instance 
in isotropic models, combine the length parameter io of edges with discrete representation 
labels as one real number, giving exponentials h{c)" = exp{inioc) = exp(mAc) with c = 
Vq^^c, X = £o/Vq^^, and /i = riA G M. In this Abelian case, homogeneous connections are 
often viewed as elements of the Bohr compactification Meohr of the real line [IQ] (rather 
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than U(l), which is obtained for fixed £o [SI])- The Bohr compactification of the real hne 

is a compact Abehan group with representations in one-to-one correspondence with those 
of M: they are given by z for all real Functions on Meohr form a Hilbert space 

using the Haar measure 



As per the Peter- Weyl theorem, all continuous functions on MBohr can be written as count- 
able superpositions 



states exp{ific) forming an orthonormal basis for all real /i. These are eigenstates of the 
derivative operator p = —Siri'-yipd/dc with eigenvalues Snjipfi. (See ST] for more 
details on the Bohr compactification of the real line.) 

The form of states suggests a map between the spaces of functions on Abelian homoge- 
neous connections and functions on the Bohr compactification of the real line: B: px^n{g) ^ 
exp{iXnc) onto the subspace spanned by all exp(i/ic) with rational fi. With the formulas for 
inner products, f[T^ and flTB]) . it follows that this map is an isometry, and it is a *-algebra 
morphism and commutes with the action of p. However, it is not invertible, and therefore 
not unitary: one can easily find (Ai,ni) 7^ (A2,n2) such that Ai^i = A2^2@ 

Not all features of the Bohr compactification are realized in homogeneous models even 
of Abelian type; care is therefore required if only the Bohr compactification is used: 

1. The label /i is a degenerate version of a pair (A, n) of state parameters, playing distinct 
roles in holonomies and discrete dynamics. The degeneracy of A and tt, in is lifted 
by a direct quantization of homogeneous connections as linear maps (j): CS — )■ CG. 

2. Our new quantization of homogeneous connections easily applies to non- Abelian mod- 
els, while the Bohr compactification of does not properly display non- Abelian fea- 
tures of general anisotropic models. Via the spectrum of our C*-algebra, we obtain 
a compactification of the space of non-Abelian connections unrelated to the Bohr 
compactification. 

The Bohr compactifiction was introduced to loop quantum cosmology in [40J by way of 
a pure minisuperspace quantization of the isotropic connection component c. Compared 

^Starting from U(l) instead of K, we make the family of representations continuous by enlarging the 
group manifold while keeping it compact. This procedure has no analog for the non-Abelian SU(2), in 
which non-trivial Lie brackets determine the representations and discrete spectra of its generators, as 
well-known from angular momentum in quantum mechanics. 

^There is a bijection between suitable subspaces of the Bohr Hilbert space and the Abelian homogeneous 
Hilbert space with structure group G = U(l). If we take the subspace restricted by < A < 1, the map 
P\,n{g) Im) |A + n) is a one-to-one transformation to the subspace of the Bohr Hilbert space with 
rational fi. With the restriction on A, one can, given /i, uniquely determine n as the integer part of fi and 
X as ^ — n. Moreover, if Ai -I- ni 7^ A2 + n2, Ai 7^ A2 or ni 7^ 712. However, choices (Ai, ni, A2, ^2) exist for 
which Ai + ni 7^ A2 + n2 but Ai^i = A27T.2; the inner product is therefore not preserved and the map is 
not unitary. 




(16) 




(17) 
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to using a periodification of the real line to U(l), as originally done in [39], this procedure 
has the advantage of faithfully representing all values of the connection component: c can 
be computed if exponentials exp{ific) are known for all real fi (irreducible representations 
of MBohr), while knowing exp(mc) with integer n (irreducible representations of U(l)), 
allows one to compute c only up to adding integer multiples of 27i. Still, this alteration of 
the original quantization is inadequate, as shown here. An isotropic connection is not a 
number c, and a diagonal homogeneous connection is not a triple of numbers c/, just as an 
inhomogeneous connection is not a collection of scalar fields. A homogeneous connection 
is a linear map from CS to CG, or an element of CS* x CG. The factor of CS* is 
crucial to relate the nature of a connection as a 1-form, but it is overlooked if one takes 
only the components C/, or a single c for isotropic models. The new quantization of 
homogeneous models provided here takes into account the correct mathematical structure 
of homogeneous connections, leading to inequivalent Hilbert-space representations. In some 
of the following sections, we will see that these differences are crucial for realizing a relation 
to the full theory and for some dynamical aspects. 

3.4 Minisuperspace operators and averaging 

Minisuperspace quantizations allow a large set of choices regarding quantum representa- 
tions, kinematical operators, and, most of all, the dynamics. The dynamics is the most 
difficult to derive from the full theory, requiring detailed projection maps to ensure that 
one stays on the space of homogeneous states; no strict derivation is available as of now. 
Fortunately, however, quantum geometry implies several general effects in the dynamics, 
for instance in Hamiltonian constraint operators of loop quantum cosmology, deviating 
from classical expressions by the use of holonomies and inverse-triad operators. The form 
of holonomies and inverse triads, in turn, is dictated by properties of the kinematical 
quantum representation used. If one can derive the simpler setting of kinematical repre- 
sentations and basic operators, properties that imply characteristic dynamics in the full 
theory are realized in reduced models as well. Reliable qualitative effects can be predicted 
even if the dynamics is not directly derived but rather constructed by analogy with the 
full theory, using reduced operators. Given that the full dynamics so far appears to be 
ambiguous, too, only generic effects are reliable, anyway. Details of the reduction of dy- 
namics may not matter much, provided one is asking the right questions. Relating models 
to the full theory helps one decide which questions can (and should) be asked. 

3.4.1 Lattice subalgebras and spin-lattice states 

For any fixed triple of integers Afj, the operators Pki/Mi,ji{.9i)^l i for all integer fc/, together 
with form a subalgebra of the homogeneous holonomy-fiux algebra, which we call a lat- 
tice subalgebra or, more specifically, the (1/A/i, l/A/'2, l/A/3)-lattice subalgebra. Any state 
Pkj/Mi,ji{.9i)^j\^)i obtained by acting with a homogeneous lattice-subalgebra holonomy on 
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the cychc state |0) independent of gi, can be written as a superposition 

Pki/J\fi,ji{9i)ni = Pi/MiJii9i)Z' P^/^fiJi(9i)tl - ■ ■ Pi/Afidii9i)m'~^ 

of products of elementary excitations pi/Mi,ji{.9i)^' ■ It can be viewed as the evaluation 
of a lattice-based spin-network state in a homogeneous connection — a cylindrical state 
whose graph is a lattice with straight edges and regular spacings = Lj/Mi- 

In order to make contact with inhomogeneous states, we use a spatial lattice of the 
form just introduced, with uniform spacing in direction Xf as measured in coordinates, 
from links along the three invariant vector fields Xf of a Bianchi I model. We require 
the region V of coordinate size Vq = L1L2-L3 to be sufficiently large, to allow many lattice 
links of the chosen spacings. We must restrict attention to Bianchi I at this stage to obtain 
closed lattices; below, in Sec. 13.5.11 we comment on other Bianchi models. For non-Abelian 
symmetry groups, such as those of Bianchi models other than type I, different generators 
do not form closed square loops by their integral curves, and therefore no lattice can be 
constructed!^ 

Fixing an orientation for each of the three directions, we label lattice links by pairs 
(f,/) of a vertex v as the starting point of a link e^j in direction (as in [S3]). For a 
connection (not assumed homogeneous at this stage), each link gives rise to a holon- 
omy hyj = V exp{J^ ^ A^Xje^ds). We will work with spin-network states of the underlying 

lattice (not required to be gauge invariant), or spin-lattice states. Each link holonomy 
appears in some irreducible SU(2)-representation with spin j^j. In the matrix representa- 
tion pj^j{h^j), we pick matrix elements {mvj\pj^j{hvj)\nyj), with two eigenstates \myj) 
and \nyj) of Pj^jir-^) (or any other component). The function {m^j\pj^ j{h^j)\n^j) is 
then an eigenstate with eigenvalues m^j and n^j, respectively, of the 3-components of 

-"^"For two generators Xi and X2 , a single closed loop is obtained if one uses integral curves of the left- 
invariant vector field of Xi and the right- invariant vector field of X2, as proposed in |48j : left- invariant 
vector fields commute with right-invariant ones. However, no complete lattice can be formed from these 
integral curves in three spatial dimensions: To generate lattice sites, one would have to fix one type of 
vector field, left- or right-invariant, for each spatial direction. If Xf is taken as left-invariant, X2 must 
be right-invariant for a closed 2-dimensional lattice in the 1 — 2-surface. For a closed lattice in the 1 — 3- 
surface, also Xg would have to be right-invariant, but then, with both and Xf right-invariant, there is 
no closed lattice in the 2 — 3-surface — unless Xf and Xf happen to commute. Lattice constructions based 
on the interplay of left- and right-invariant vector fields cannot be performed for all Bianchi types, making 
those constructions in the available cases (Bianchi I and II) non-generic. Attempts at such constructions in 
anisotropic models show some of the pitfalls of ad-hoc assumptions, as illustrated by the series [49 1 150 1 151] 
of papers where most initial claims of (49) . for instance regarding averaging or a possible relation to 
lattice constructions, had to be withdrawn or weakened in later installments. Initially simple-looking 
constructions became more and more contrived. Instead, it is more general to use lattices according to 
the kinematical structure of Bianchi I, and then implement other Bianchi models by suitable curvature 
terms in the dynamics [52]. In this way, all Bianchi class A models can be quantized with one and the 
same scheme. One may worry about an inconsistency in using Bianchi-I lattices for other Bianchi models. 
However, at the inhomogeneous lattice level, no strict Bianchi models can be realized. The symmetry type 
just provides guidelines along the way to consistent dynamics, which can well be realized for all class-A 
Bianchi models. 
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right-invariant and left-invariant derivatives by h^j. Our spin-lattice states are therefore 
functions V'(j„,7,m„,j,n„j)(/i) = Ilyj{mv,i\Pj,j{K,i)\nv,i) depending on the connection via 
link holonomies, with an inner product defined as usual by integrating over all h^j using 
the Haar measure |5l]. This defines the Hilbert space T^iattice- For short, we will write 
these states as \{jv,i,'m'v,i,nvj)). We have the usual action of holonomies and fluxes. 

3.4.2 Homogeneous distributions 

A homogeneous analog of spin- lattice states, depending on holonomies h^^XjTi) = exp(A/0(T/)), 
is ip{\i,ji,mi,ni){gj) = Ili{mi\pj{h^{XiTi))\ni) = Yli{mi\pxjjj{gi)\ni) with A/A/} integer, 
written for short as |(A/, j/, mj, nj)). There is an additional label A/, replacing the edge or 
link dependence of inhomogeneous states and representing the CS-pait X = XjTj of a ho- 
mogeneous connection G £S'*®£SU(2), subject to the scaling condition. The set of these 
states is fixed by holonomies in the (l/Ai, 1/J^2, l/A/'3)-lattice subalgebra, with elementary 
holonomies acting by multiplication, changing the SU(2)-representations ji according to 
recoupling rules, and flux operators having eigenvalues STT'yipmj (for right-invariant vector 
fields) and STT'jipni (for left-invariant ones). No decomposition as in (|6]) is required since 
we have a fixed common denominator A/} for all holonomies considered in direction Xj. 

So far, homogeneous and inhomogeneous lattice states are defined separately from 
each other. We relate them by introducing a map a: "Hhom — ^ lattice, \ {Xi, ji,mi,ni)) i— )► 
((■^Z) J/) \ from the homogeneous Hilbert space to distributions on the lattice Hilbert 

space. This map is the key ingredient of quantum symmetry reduction, as described in 
Section [21 Following [19], we define homogeneous distributions by their evaluations 

((A/, ji, mi, m) I (j„,/, m^j, n^j)) = ((A/, j/, m/, m) \ {j^j, m^j, n^,/)) U„ ,=exp(0(r,)/Ar,) (18) 

on all basis states of ^lattice- On the right-hand side, the inner product is taken in T/hom, 
with \{j ,i))\Kj=exp{4>(Ti)/Ni) obtained by restricting the connection dependence 

of the spin-lattice state to homogeneous 0. The distributional evaluation vanishes unless 
the representation ji appears in the tensor product ^^jjvj, and m/ = Ylivi''^v,h "^i = 
The reduction of states depends on the size of the region V via A/}, just like 
the classical reduction of the phase space. 

3.4.3 Averaged operators 

An operator O can be reduced from the lattice theory to the homogeneous Hilbert space if 
its dual action fixes the space of homogeneous distributions: If there is a l^/^) G Hhom such 
that {{Xi,jj,mi,ni)\d^{j^j,myj,n^j)) = {il)\{jv,i,m^,i,n^,i)) for all Kj^,,/, m^,,/, n^,/)), we 
define 0|(A/, j/, m/, n/)) = All link holonomies {rriy^jlpj^ j{h^ j)\ny^j) along symmetry 

^^In (fT8|) . we restrict to holonomies hy,i = exp(A/(/)(T/)) with u-independent A/ = 1/Afi, or a regular 
aligned lattice of uniform link lengths. At this stage, we could allow irregular lattices with varying £i{v), as 
long as all links are still along symmetry generators XiTi. Different lattice sectors would then contribute 
to the reduction, and refinement would be necessary in the multiplication and action of holonomies. This 
option will be discussed in more detail below. 
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generators, taken as multiplication operators, satisfy this condition. They act on distribu- 
tional homogeneous states by {Tn^j\pj^ j{K^j) {gi\i\i, ji, mi, m)) = pi/Afj,j,,jigj)n^,j ■ 
PXijiidi)^' , just as in the reduced space of homogeneous states. 

Flux operators require additional constructions. A single lattice flux F^j associated 
with a surface dual to link Cyj does not map a homogeneous distribution to another such 
state: Take a set of states li'vj) '■= |0, . . . , 0, (1/2, 1/2, 1/2), 0, . . . , 0), each with non-zero 
labels only on one lattice link Cyj. We have 

{{l/Ml,jl,mi,ni)\Fyj\lljyj) = 47r745ij,l/2^mj,l/25nj,l/2 

and {{Xi, ji,mf,nf)\Fyj\ipy'j/) = if f 7^ or / 7^ /'. Therefore, 

((1/M, 1/2, 1/2, l/2)\Fyj\iljyj) ^ ((1/A/}, 1/2, 1/2, 1/2) IF^,,!^'.',/) for v v' . 

However, we must have {"^lipvj) = {"^li^v',!) for any homogeneous state 6 "Hhom 
since V'»;,/U»,/=exp(Aj?i(Tj)) = V'i>',/U„,i=exp(Aj^(Tj))- The state ((1/M, 1/2, 1/2, 1/2)| cannot 
be contained in a decomposition of {{1/J\fj, ji,mj,nj)\Fvj in our basis, and we can re- 
peat the arguments with arbitrary values of the non-zero label in lipvj) to show that no 
homogeneous state can be contained in the decomposition. Therefore, the distribution 
{{l/Afi,ji,mi,ni)\Fyj cannot be a superposition of homogeneous distributional states: 
flux operators associated with a single link do not map the space of homogeneous states 
to itself. 

Even classically, the flow {■, Fs} generated by a flux operator is not everywhere tangent 
to the submanifold of homogeneous connections and unrestricted triads in the inhomoge- 
neous phase space, but it is tangent to the subspace on which both the connection and 
the triad are homogeneous. In the quantized theory, using distributional states, we have 
ensured states to be restricted to homogeneous connections, but no such condition has yet 
been implemented for the densitized triad or fluxes. 

Flux operators must be averaged to generate a flow that keeps the space of homogeneous 
states invariant. However, non-Abelian gauge properties prevent us from simply adding 
J2n Is ^i^ad'^y for a family of surfaces Sn translated along the generators of the symmetry 
group. Instead, we must relate the fibers of the SU(2)-bundle in which takes values, 
using parallel transport between the S'„. (This problem seems to be related to issues 
encountered in constructions of a non-Abelian triad representation [55] . Here, homogeneity 
will help us to propose a solution.) 

To describe the specific construction, we assume an aligned state, consisting only of 
holonomies h^j in the three independent directions but not necessarily forming a regular 
lattice. For an averaged , we choose families of surfaces Sn,i transversal to the symmetry 
generators Xf, such that they have co-normals = 5^ , layered at regular intervals across 
the region V. Eventually, we will send the number N of surfaces to infinity. Before doing so, 
we define a gauge-covariant averaging by = A^""*^ 'l2n=i Is i ^dhi{y){Ei{y)nl)(}?y where 
hiiy) is the connection-dependent parallel transport from some base point, chosen for each 
integral curve in direction J, to a point y on the surface. 
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The base points will be chosen in a state-dependent way because the state determines 
how the connection is excited, usually in a discontinuous way at lattice vertices. We 
decompose a state as a superposition of contributions ^ = 'ipih^^j, h.uji^K)Ylyj Pj^. i{hvi,i) 
where the dependence on holonomies along directions J and K will not matter. The 
set of all vi then gives us all vertices where parallel transport in the /-direction changes 
discontinuously. We will average with these vertices chosen as base-points, so that only the 
continuous parts of parallel transport are taken into account. We first decompose surfaces 
Sn,i = IJfc Is J fc thai each piece Sn,i,i intersects at most one edge. In the action of the 
fiux operator, instead of summing over k we will then be summing over edges intersecting 
the surface: We write 

1 ^ r 1 ^ 

n=l ''^^,1 n=l e„^,jn5„,j^0 

with h^jj(vn) the parallel transport from vj to the intersection point Vn of e^^j with Sn,i- 
If a piece of the surface Sn,i intersects an edge e^jj, we chose the base point to be 
vi (for a right- invariant vector field quantizing the fiux, or the other endpoint of e^jj for 
a left-invariant one). The adjoint action of h^jjivn) in the averaged fiux then implies 
that a fiux operator does not insert just Tj in the holonomy, at the intersection point 
{vn} = Snj^ie^jj with the surface, but hyjj{vn)^^Tihyjj{vn)- For a single edge e, splitting 
the holonomy he := h^\vn)h^\vn) in two pieces h^\vn) and h^\vn) at an intersection 
point Vn, we thus have 

1 ^ 

pl{m'\pj{he)\n') = -STTijil lirn^ — ^(m'|pj(/i(^)(t;„)ad^(i)(^^^(Ti)/ii^H^^n))|n') 

n=l 

1 ^ 

n=l 

where only terms such that Snj fl e 7^ contribute. For large A^, the number of non-zero 
contributions divided by N approaches Ag, the ratio of the length spanned by he relative 
to Lj. We obtain eigenvalues STrjipXein' . Not surprisingly, in this homogeneous context 
the averaged fiux does not refer to any point on the edge where it acts, but it picks up the 
relative length of the edge by the number of intersection points. For multiple edges, the 
fiux acts by the product rule. 

If all edges involved form a regular lattice, with the number J\fj = Lj/ij of lattice links, 
it follows that 

^' = :y^E^V.- (20) 

The factor of l/A/} = ii/Lj eliminates over-counting by adding fiuxes of all lattice sites 
along direction /. In the other two directions, on the other hand, we sum rather than av- 
erage because the minisuperspace p[ is defined for a surface stretching through the region 
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V, as in ([8]). Indeed, heuristically, the eigenvalues of pj, Sn'yipAfj ^ J^v^^J '^^^ writ- 
ten as LjLk multiplying the average value of the densitized triad: 8n'y£pJ\ff^ J2v ''^vj = 
STTjipAfjMKmJ = AfjAfxE^, where mj = (A/iA/'2A/3)~^ m^,/ is the lattice average, 
quantizing the average of the plaquette flux J E-^cix'^dx^ /STi^yip. 

3.4.4 Quantization commutes with symmetry reduction 

Holonomy operators in T^hom are directly obtained from their dual action on ©lattice- It 
does not matter whether we act with a holonomy operator first and then symmetry-reduce, 
or first reduce and then act with the corresponding homogeneous operator: 

a{pi/Mjj,^ j{gj)\{Xi,ji,mi,ni))) = pj^ j{Kj)a{\{Xj,ji,mi,ni)) . (21) 

After averaging, the same commutation relationship is realized for lattice flux operators. 
We have 

ii^/^i,ji,mi,ni)\fi\{jyj,myj,nyj)) = ^((l/TV}, j7,m/,n7)|F„,j,3|(j,,,7,m„,/,n^,/)) 

87r74 ^ 

~ / ,''^v,J ^ji(^®„jv,i^rni,Y,„m^,i'^ni,Y.„'n-v.l 

= 87174 Ajmj i{Xi,ji,mj,ni)\{jyj,m^j,n^j)) 

using the ^-identifications and Aj = 1/A/j in the last step. On reduced states, on the 
other hand, we have p^\{Xi, ji,mj,ni)) = Sn'-fipXjmj from the right-invariant vector field 
dH]) in the Aj-sector. Comparing these equations, we see that a{p'l\{Xi, jj,mj,ni))) = 
p^a{\{Xi, ji,mi,nj))), with analogous calculations for other components of pf . (As re- 
marked in Sec. I3.3.2[ the non-Abelian pf is not a symmetric operator unless we are in the 
lattice setting of fixed A/ for all states involved. Similarly, p/ is not symmetric in this 
situation due to averaging, in particular with a state-dependent choice of base points for 
parallel transport.) 

For basic operators, it does not matter whether we quantize or reduce first. We obtain 
the same representation properties as in the full theory, and the same qualitative quantum- 
geometry effects. But a quantitative correspondence is more complicated for composite 
operators, especially the Hamiltonian constraint crucial for dynamics. 

3.4.5 Holonomy-fiux algebra in reduced models 

The previous calculations have shown how holonomy-fiux representations of homogeneous 
models are derived from the full algebra. Using the minisuperspace embedding a, we 
obtain basic operators — holonomies and fluxes — from their action on inhomogeneous 
lattice states. Since the holonomy-fiux representation of the full theory is unique [561 
IFF] , the minisuperspace representation derived here, by restriction of the full algebra to 
lattices followed by taking the dual action on homogeneous distributions, enjoys the same 
distinction. 
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So far, we have written all formulas for the general case of non-Abelian homogeneous 
models. Using the same techniques of restriction of states and reduction of operators, it is 
straightforward to implement diagonalization or isotropy: We restrict states to (^} = C(/)A} 
for diagonal models, or 0} = cA} for isotropic ones. Flux operators = A'^^pf then leave 

diagonal states invariant (while Ajpf with J ^ I would not), and the averaged V = \ J2iP^ 
leaves isotropic states invariant. These situations are covered in |35j . 

3.5 Dynamics 

From holonomy and flux operators, we construct more complicated ones such as the volume 
or the Hamiltonian constraint. The volume on spin-lattice states can be defined as in 
the full theory, using \\e^^'^e^'^^ F^j^iF^^jjFy^K,k\^^'^ ■, just restricted to 6-valent vertices as 
encountered in a lattice. The complete spectrum is unknown in the non-Abelian case. For 
simpler algebraic relations, we may replace the cubic SU(2)-invariant with a product of 
quadratic invariants, K := Y\^i=i{Fv,i,iFv/Y^'^ with eigenvalues (87r7)^/^£| Y\j=i{jv,i{3v,i + 
1))^/^. In what follows, details and differences of these spectra will not play a major role, 
and we will make use of the simpler version. 

As part of the formulation of dynamics, we will be interested in reducing the vol- 
ume operator. An important question for non-linear combinations of basic operators 
is whether the average is taken before or after reduction. The minisuperspace volume 

I{]=Mf''f'^ = {^1^2^:,)-^'^ nLi using Jj = and (EO]), has eigenvalues 

n </J^u^) = (8-7)^/^4 n v^imiTTT) . 

The spectrum can be computed using pj in (jS]) or p[ in ( l20l) . but it does not equal that 
of the averaged {J^iJ^2-^3)~^ "^yVv, which would be the reduced volume operator. The 
distinction is important for the correct size of quantum-geometry effects, as we will see 
below. Pure minisuperspace models make use of pf or p[ ; correctly capturing quantum 
effects requires an averaged volume operator. 

3.5.1 Hamiltonian constraint 

The classical Hamiltonian constraint contains curvature components, to be represented 
in homogeneous models by holonomies P\,j{gi)^, non-polynomial functions of the con- 
nection which differ from connection or curvature components by higher-order terms. In 
Abelian models, it is is easy to see that px^nidi) as an operator is not weakly continuous 
in A at A = 0, and it is not possible to define a connection operator via d/dA|A=o: The 
diagonal matrix elements {{X'j,nj)\px^n{gi)\{^'i,n'j)) = ^o.An, using ([T^ . are not continu- 
ous in A at A = 0. In the non-Abelian case, the argument is more complicated because 
{{X'j, jj,m'j,n'j)\pxj{gj)^\{X'j, jj,m'j,n'j)) need not be zero for A 7^ 0. The value rather 
depends on the multiplicity of the trivial representation in the tensor product pf,'^'^ ® pf^ 
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if = q/r and \ = p/r with their least common denominator r. For A = 0, we pick out 
the trivial representation in p^^''; for A 7^ we pick out the coefficients of all irreducible 
representations of p®^ in p^^^- In general, these coefficients are not continuously related 
at A = 0. 

Instead of using derivatives and connection operators, we are required to use holonomies 
PsjiOi)^ with some finite 6 to construct the Hamiltonian constraint operator, appealing 
to the standard relationship between curvature components and holonomies around small 
closed loops. Since the same relationship is used for the full constraint [581 ES], there is at 
least a plausible connection between models and the full theory. We have 

hA = l + fsls'AT, + Oit) (22) 

if the loop A, spanned by two unit vectors and S2, is of coordinate area i"^. In lattice 
constructions, one may use loops around elementary plaquettes, although consistency is- 
sues of the constraint algebra may require more complicated routings (see e.g. [33]). The 
coordinate area of a loop in the (/, J)-plane is then close to iiij, and we are led to use 
homogeneous holonomies psjjjigi) with 61 = ii/Lj = 1/A/}. 

In Bianchi models, we have the Hamiltonian constraint ([2]), or ([3]) when diagonalized. 
Putting holonomies along square loops and suitable constructions of triad operators to- 
gether as in the full theory [59], we obtain, following [371 138], 

^ " ~ (87r)^^^G^h6 6 6 ^ e"'^t'^(p5/,i/2(^7/)P5„i/2('7j)P5„i/2(^/)"V<5^,i/2(^j)'' 

X P5j„i/2igK)[p5j„i/2igKy\V]\) + Hr (23) 

with matrix products of all holonomies involved. Instead of j = 1/2 one may use other 
irreducible representations, or add different such contributions [6OI [611 [62] . Also the values 
of 61 are subject to choices whose implications we will discuss in more detail below. One 
may choose fixed values, or relate them to properties of the state acted on by the constraint 
operator. In the latter case, the state dependence of regularized constraints in the full 
theory would be modeled. For now, however, we will assume fixed 61 so that H preserves 
a lattice subalgebra. 

The term Hr in f l23]) vanishes for the Bianchi I model and incorporates spin-connection 
terms for other models, as in [52j. The commutator \psj^^i/2{gK)[p5K,i/2{gK)^^,V]\ in the 

second line quantizes the classical combination e^^^E^E^^ I \/ \ det E\ which diverges for 
degenerate triads, at classical singularities of collapse typeo Writing this operator on 



^^The absolute value around the commutator is necessary because the classical analog = 
27r7Gsgn(det Ey Ef E^^ / ^\AGi E\ carries a sign factor. The absolute value avoids parity violation 
(see [53] for a detailed discussion of parity) . Classically, the sign of det E changes whenever the sign of 
{A^,V^} changes, and one could multiply the latter with sgn(deti?) to avoid parity violation. However, 

when quantized, P5^,i/2(<?i<')[P(5^.i/2(5K)~^, ^] and det do not commute in non-Abelian models. Since 
the dynamics identifies states annihilated by the commutator as degenerate ones corresponding to classical 
singularities, we refer to its own sign instead of multiplying it with the non-commuting operator sgndet E 

m. 
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the right in H, in the ordering chosen in (123!) as in the full theory [59], implies that 
singular states decouple from the dynamics: They are automatically annihilated by the 
constraint. (There is a similar factor in ifr, also ordered to the right.) In this way, quantum 
hyperbolicity [7] is realized even if no difference equation is available as an explicit evolution 
equation. 



3.5.2 Abelian models and difference equations 

Difference equations are obtained in Abelian models after transforming to the triad rep- 
resentation, provided the Hamiltonian constraint fixes a lattice subalgebra. First, writing 
Abelian holonomies of a diagonal Bianchi model as psj,i/2{gi) = cos(5/C//2)+2A}rj sm{SiCi/2), 
one can compute all matrix products and the trace in ( !23l) . Since the final result is lengthy, 
we define 

K3 := sin((5iCi/2) cos((5iCi/2) sin((52C2/2) cos((52C2/2) = - sin(5iCi) sin((52C2) 
and cyclic permutations thereof, as well as 



2i {sm{6jCj/2)Vcos{6jCj/2)-cos{6jCj/2)Vsm{6jCj/2) 



The Hamiltonian constraint is then 

H 



1 



(24) 



j=i 



It is straightforward to compute the action of Kj and the eigenvalues of Ij on the (5i, 52, ^3)- 
lattice subalgebra (now dropping the fixed 5/ from the notation of states): 



k^\ni,n2,n^) 



-777(1^1 + 25i, n2 + 2^2, r^s) - |^i - 25i, ^2 + 2^2, ^3) 
-\ni + 25i,ra2 - 2^2,^3) + \ni - 26i,n2 - 262, ris)) 



(25) 



and cyclic permutations, and 



/i|ni,n2,n3) 



I ^^i+<5i ,n2,n3 ^ni—Si ,n2,"3 11^1' ^^2 j ^3) 



(8717)^/24 ^\n^ + s^\- ^\n^-s^\ v/|n2n3||ni,n2,n3). (26) 



Since H\ni) must vanish (or equal the action of a matter Hamiltonian -^matter on |n/)), 
a difference equation is obtained in the triad representation of coefficients ipni,n2,n3 = 
Y.ni^n^,n2,n3\ni). Introducing s„,,„2,„3 := VlnI^i2^V'ni,n2,n3 to shorten the expression@ 



^■^Dividing by nin2n^ is well-defined for the evolution equation because |ni, 712,713) is annihilated by 
the constraint whenever 711712713 = 0, as part of the property of quantum hyperbolicity. The coefficients 
'4'ni,n2,na with 7ii7i27i3 = dccouplc from the rest and can safely be ignored. 
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we have 



^Sii^l) (-^m, 712+2(52, n3+2(53 ■^^-^^^j— 252,n3+2<53 ^ni,n2+252, 713—283 ^ni,n2—252,n3—2S3) 

^52(^2) ('5ni+2(5i,n2,n3+2<53 ^ni—2Si,n2,n3+253 '5ni+25i,n2,n3— 253 ~^ '^«i~2(5i,7i2,n3— 253) 

^53(^3) ('5ni+2(5i ,712 +2(52, n3 "5711— 2(5i,n2+2(52, 713 '5ni+25i,n2— 2(52, 713 ~^ '^ni—2(5i,rj2— 2(52, 713) 
10Q^2„3^/)2 r r r -^matter (""-l, ""-2, ^3) ^ 

1287r 7 G£pdid2d3 z- Sni,7i2,n3 (27) 

^7ii,n2,"3 



with volume eigenvalues inserted. We used partial eigenvalues of the matter Hamiltoniaio 

-^matter I ^) = ^ (^-f^matter (^/) 
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which may still act on a matter-field dependence of s„j „,2,n3- The other coefficients refer to 
eigenvalues of //, defining Asj{nj) := Ii{ni,n2,n3)/Vn-^,n2,n3 such that As{n) = \\/\n+~5\~ 
^/\n^^^\/ ^y\n\. Equation (127|) correctly quantizes the classical terms such as 020301 = 
C2C3^y\j^^Pjp^ in the Hamiltonian constraint ([3]), using the ordering | ""^0203 a/ \p^p'^p^ \ . 

Equations such as ( !27|l have been derived long ago [38], and reproduced since then 
many times, in slightly different forms. We have presented the derivation here with some 
detail because, lacking the proper notion of homogeneous connections, it had not been 
realized before that the form is valid only for a Hamiltonian constraint operator fixing 
a lattice subalgebra of the homogeneous model. If the 61 are not fixed once and for all 
but depend on the state acted on, or if the 5/ used in the operator are not the same as 
those of the lattice subalgebra, the decomposition rule (EI) or its Abelian analog ( 1T3|) must 
be used, implying refinement. The operator ( 123|) will still be valid, but its action must 
be re-derived, and does not easily give rise to a difference equation, certainly not one of 
constant step-size. We emphasize that refinement is realized even if the operator ( 123|) is 
not modified, provided only one applies it to states not in the lattice subalgebra fixed by 
it. 



3.5.3 Lattice refinement: a toy model 

So far, we have presented the minisuperspace quantization of Hamiltonians. For contact 
with the full theory, we must try to reduce an inhomogeneous operator and face the aver- 
aging problem. This task, at present, cannot be done in detail, but its outcome will affect 
the choice of 61. Instead of deriving these values and their relation to states, we must 
resort to sufficiently general parameterizations to model different possible reductions. 

If the 61 are not adapted to a lattice or to the common denominator of all holonomies 
involved, H will not fix any lattice subalgebra, even if the 61 are fixed and not state- 
dependent. Lattice refinement then occurs by multiplying holonomies of different edge 
lengths and obeying the decomposition rule Evaluations of the Hamiltonian constraint, 

^^For simplicity, we assume the absence of connection couplings; see [65] for a fermionic model in which 
the assumption does not hold. No qualitative changes to the present statements occur in such a case. 
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especially the non-Abelian version, become more involved and difference equations no 
longer are readily available, but the property of lattice refinement can be seen already in 
an admittedly rough toy model. 

Let us assume that each 5/ is always half the maximum A/ encountered in a lattice- 
subalgebra state acted on. In terms of an inhomogeneous lattice, this means that every 
new edge generated by a vertex contribution of the Hamiltonian constraint goes half-way 
to the next vertex. In an inhomogeneous lattice, the presentation of refinement depends 
on the order in which individual holonomies or vertex contributions of the Hamiltonian 
constraint act. If refinement proceeds regularly, staying close to cubic lattices of nearly 
constant link lengths, one would expect that all plaquettes will first be split half-way along 
edges, and when this has happened for all of them, one would proceed to quarters and so on. 
However, in the Hamiltonian constraint all vertex contributions appear in superposition, 
not in simultaneous action on a single lattice. To realize an ordering, one may assume that 
a physical state annihilated by the constraint is expanded in spin-lattice states according 
to the eigenspace of some operator such as the total volume, by the maximum spin on all 
edges, or by the number of plaquettes. Ordering spin-lattice contributions in a physical 
state with respect to any of these values, plaquettes will be filled in a certain arrangement, 
such as the one described. 

Back in our homogeneous model, starting with a state in some lattice subalgebra with 
values xf\ the first action of the Hamiltonian constraint, multiplying with holonomies of 
lengths Xf^ /2 in different directions, requires a decomposition ([6]) of the whole state, refin- 
ing edge lengths to halves. After a single multiplication, no homogeneous holonomy of the 
original length xf^ will occur explicitly, but such edges are still present in a corresponding 
homogeneous lattice because there are several matrix products Px'-'^^ /2 P\'-'^^ /2 j^9i)n 

of two Aj'^''/2-holonomies without intermediate factors. Taking these holonomies into ac- 
count, we keep acting with A^^V^-holonomies until all those products disappear. (In in- 
homogeneous lattice language, we fill all Aj°''-plaquettes with new edges and vertices at 
all midpoints of the original edges.) Once these options have been exhausted, the next 
refinement step is due, going to X^P /A until all the previously refined plaquettes have been 
filled. In the process just described, we have assumed a certain ordering of the actions of 
individual vertex contributions, first filling all the Aj''''-plaquettes, then moving to Ai°V2- 
plaquettes, and so on, as we would do on an inhomogeneous lattice. 

We can relate the number of plaquettes, or the degree of refinement, to geometrical 
quantities of the whole lattice. Starting with a nearly homogeneous lattice with all edge 
spins jo equal, the initial area in the J, -ft^-plane is approximately Aq = (Aj°''jo)^A/jA/i^, 
with MjMk plaquettes in this plane and transversal links of the size Xf^ = Lj/Mi. Thus, 
Aq = L'^jIMjN'k / 1 or Aq Vq^^Jq for a nearly isotropic lattice with equal edge numbers 
in the three directions. When all these plaquettes have been refined after the first stage, 
the maximum spins have changed to ji = 2jo + 1/2: multiplying jo with two because 
we decompose holonomies halfway according to ([S]), doubling them over, and adding 1/2 
from the action of a new holonomy in the fundamental representation. The added 1/2 
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will soon be irrelevant when j becomes larger by repeated doubling. The area has then 
increased to Ai = {\^pf{2iQf{2Mj){2MK) with A^^^ = /2. Combining these equations, 
Ai = 4:Ao, in which only the increased spin due to refinement contributes. After N steps, 
the same arguments show that the area has increased to An = 2'^Ao, and A^^'' = xf^ /2'^ = 
^y Aq/ An- The spin quantum numbers increase by Jn ~ 2^jo- 

It is an interesting feature that the spin of a single action (here 1/2), an ambiguity 
parameter of the full constraint, becomes progressively less important in the model as 
refinement proceeds, increasing j^. The large-scale behavior is insensitive to details of 
the microscopic dynamics and associated ambiguities, a property that makes effective and 
mean-field viewpoints meaningful. 

In this model, the edge lengths A^^^ are inversely proportional to the square root of the 
area, or to a linear measure of the extension of the lattice. With near isotropy, this scaling, 
X = Xo/a, is of advantage for holonomy-modified dynamics, in which holonomies exp(iAc) = 
exp(iAoc/a) depend on the isotropic connection component c only in the combination 
c/a (X proportional to the Hubble parameter; the same behavior has been proposed in 
[55] as an ad-hoc choice. While c may grow large even at small curvature, for instance if 
there is a positive cosmological constant, H remains small in low-curvature regimes. The 
refined dynamics, with a non-constant A, is more well-behaved in semiclassical regimes. 

In diagonal anisotropic models, the refinement behavior described here implies that 
holonomies depend on the connection by the combinations cjj ^/\^\^ as in [67], being 
proportional to the area A of the plane transversal to the /-direction. This refinement 
is problematic in terms of stability properties of the difference equation it implies [68] . 
A different refinement scheme in which A/ is inversely proportional to the length of the 
/-direction is preferable; a more advanced action of the Hamiltonian constraint not fixing 
a lattice subalgebra would be required, a refinement scheme in which 5/ depends not only 
on the length A/ of its own direction but also on the other two links meeting at a vertex. 

3.5.4 Difference equations with mean-field refinement 

Lattice refinement is the homogeneous realization of discrete dynamical processes in the 
full theory; ideally, its form would be derived by reducing a full Hamiltonian constraint. 
A dynamical state of quantum gravity should in general be expected to have different 
lattice structures and spacings at different times, or on different spatial slices, especially in 
loop quantum gravity whose Hamiltonians are generically graph-changing. The number of 
lattice sites is then a dynamical parameter. Indeed, if the A/} or the li are kept constant — 
we always assume Li, L2 and L3 to be constant as these are classical auxiliary parameters 
— cosmic expansion would quickly blow up the discreteness scale, ^i\/\p^'p^v'^\l\v^\ as 
measured in a diagonal Bianchi geometry, to macroscopic sizes. Lattice refinement must 
be a key feature of quantum-gravity dynamics. Dynamical minisuperspace operators such 
as the Hamiltonian constraint should not refer to constant A/} or 5/, as assumed so far in 
fl23|l and fl27|) . but to parameters that depend on the total volume or the scale factor via 
an evolving discrete state. 

Strict difference equations of loop quantum cosmology then do not exist, even in Abelian 
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models, and approximations cannot always be derived easily. The correct evolution equa- 
tion in a triad representation would rather have to implement the changing number of 
degrees of freedom, a problem studied in other contexts as well [Ml EQI [71]. Instead of 
working with such complicated equations, there are two approximation schemes that help 
to find properties of solutions: Effective equations and difference equations in redefined 
variables. 

Effective equations describe properties of solutions of difference equations in Abelian 
models via a non-canonical basic algebra, the discreteness implemented by using 
exponentials of the connection. For instance, we would represent a discrete degree of 
freedom {c,p) by a non-canonical basic pair {exp{i6c),p) with a closed linear algebra 
under Poisson brackets. If 6 depends on p by a power law 6{p) = So\p\^ as a form of 
lattice refinement, {exp{iSo\p\^c), \p\^~^) still satisfies a closed linear algebra [72]. We 
then generate evolution by a Hamiltonian much like f l23p . depending on exp{iS{p)c) 
according to the regularization chosen. Effective quantum evolution equations then 
follow the general scheme of [731 El] and provide approximate information about 
refining solutions. While strict difference equations are not available in lattice-refining 
Abelian or in non-Abelian models, effective equations can still be formulated and 
solved in both cases. 

Difference equations in re(de)fined variables model refined quantum evolution by 
difference equations equally spaced not in the original triad eigenvalues, but rather 
in some redefined versions obtained as non-linear function of them, such as power 
laws. One can derive a suitable equidistant parameter if one knows how the 5/ de- 
pend on rij in the refining case. Instead of an equation (1271) with rij- dependent 
increments, for instance in Sni,n2+2S2ini,n2,n3),n3+2S3(ni,n2,n3), one can sometimes work 
with an equidistant difference equation in re(de)fined independent variables. If 6i 
depends only on rij with the same value of J, we define fiiijii) := j^\5i{z))^^dz 



such that ni{ni + 5i{ni)) = j^'^''^'''\5i{z)Y^Az = n,(n,) + = 



niijii) + 5i{ni){5i{ni))^^{l + 0((5^(n/))) = niini) + 1 + 0((5^(n/)), a constant incre- 
ment in regions in which the derivative is sufficiently small. If oc 
is a power law with x < for refinement, the equidistant approximation is good 
at large n/ but not for small n/, where the quantum dynamics remains ambiguous, 
anyway. (For 5/ oc we have an equidistant equation in fij oc |n/|"^~^, corre- 
sponding to the new p-dependent variable used for effective equations.) One may 
also redefine the whole difference equation in terms of ftj with constant increments, 
dropping 0(5j(n/))-terms as a specific choice of factor ordering [68] . 

If 5i depends not just on rij with the same /, a redefinition is more complicated 
to derive. If 5i525^ is proportional to a power of \nin2n^\, 81828^ oc \nin2n^\^ such 
that refinement does not introduce additional anisotropy, one can always find one 
equidistant variable given by 




{5i{zi, Z2, Z^)52{zi, Z2, Z^)5z{zi, Z2, z^)) dzidz2dz:>, 



30 



oc \nin2n3\ 

related to the total volume |68j. This choice resembles Misner variables [75], which 
refer to the volume (or scale factor) and two anisotropy parameters. 

A state dependence of dynamical operators, underlying lattice refinement, may seem 
unexpected. After all, the dynamical operators are used to derive evolving states; how 
can properties of such states enter the definition of dynamical operators or the difference 
equations they imply? Taking reduction seriously, it turns out that state dependence is 
unavoidable. In minisuperspace models, we cannot formulate a dynamical operator from 
first principles, or if we do so, the results are fraught with minisuperspace artefacts because 
full properties of the discreteness are ignored. As described before, reduced dynamics is 
supposed to model the full dynamics of a symmetric state, to be projected back to the space 
of symmetric states after each application of the evolution operator, a process that includes 
the decomposition rule ([6]) used crucially in our toy model of refinement. A minisuperspace 
evolution operator obtained by reduction must encode both the full Hamiltonian constraint 
and properties of the projection. The latter depends on the evolving state to be projected 
back on the symmetric space. While the precise form remains complicated to determine, 
we see how a state dependence of the end result is obtained. 

Without a detailed method to perform dynamical reduction, the phase-space depen- 
dence of parameters such as 6j or A/} is inserted in equations only after the operator has 
been formulated, as a kind of mean field describing microscopic properties not directly ac- 
cessible at the minisuperspace level. As already noticed, large-scale dynamics is insensitive 
to microscopic details. Such details and ambiguities are relevant at small scales and at 
higher curvature, or in strong quantum regimes. These regimes can be understood only 
by general effects, such as quantum hyperbolicity, but away from deep quantum regimes, 
effective and mean-field pictures are meaningful and useful. 

3.5.5 Ad-hoc modifications 

We have presented a quantization of homogeneous models which has a tight link with the 
full theory and, unlike previously existing versions, applies in non-Abelian cases. Lattice 
refinement naturally arises as a consequence of state-dependent regularizations as in the 
full theory, combined with a reduction of all states, including physical ones, to the space 
of homogeneous loop quantum cosmology. 

Lattice refinement is important for consistent dynamics, for a fixed lattice expanded 
by cosmic evolution would either be coarse at the present time, or would have to start at 
tiny spacings, orders of magnitude below the Planck length, to be unnoticeable in current 
observations. Lattice refinement as a dynamical process ensures that the discreteness 
scale does not need to follow cosmic expansion; it can remain small at a constant or 
slowly- changing value as macroscopic events happen on larger regions. With a concrete 
realization of lattice refinement, we can look back at minisuperspace modifications that 
have been proposed in the hope of obtaining appropriate dynamics, and see how justified 
their assumptions are from the perspective of the new picture. The most commonly used 
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ad-hoc modification is a change of classical basic variables before isotropic minisuperspace 
quantization, in which c/a oc a/ a appears in holonomies, and the role of p is played by the 
volume. 

We first note that fluxes necessarily result as reduced operators in the derived basic 
algebra, not other powers of densitized-triad components or the volumel^ Basic operators 
or linear functions of them are directly reduced by reference to the commutation result 
of Sec. 13.4.41 Non-linear functions such as the volume, on the other hand, are more 
complicated to average or reduce exactly, with no currently known procedure to do so0 
The basic holonomy operators therefore act by shifts on the flux spectrum by constant 
amounts (of p in isotropic models), not the volume spectrum. If constant shifts of the 
volume spectrum have dynamical advantages, as in the model of [66], they cannot be 
derived by direct use of basic operators but only after re(de)fining variables as in the 
preceding subsection. The volume can be used as a basic variable only as a modification 
within a pure minisuperspace quantization, without reduction and a justified analog in the 
full theory. 

The modifications proposed in [66] have been motivated in holonomy-based expressions 
for in the constraint by referring to geometrical areas a^^'^ instead of coordinate ones 

where a is the scale factor of a Friedmann-Lemaitre-Robertson- Walker model to be 
quantized. However, just as tensor components depend on coordinates, it is the coor- 
dinate area £^ which should be used in the expansion ( l22l) . not geometrical areas obtained 
using the metric or densitized triad. (Contracting with the two vector fields provides 
a scalar. However, for the coordinate area £^ to be the correct factor in the expansion, the 
vector fields must be normalized using a background metric. Changing coordinates and 
retaining normalization then makes the contracted version transform.) 

If these and other ad-hoc assumptions, for instance about factor ordering, are dropped, 
dynamical equations are much more ambiguous than usually realized or admitted. More- 
involved constructions of lattice refinement are required, which capture necessary projec- 
tions of the dynamical flow back on the space of symmetric states. Exact projections being 
largely unknown, the dynamics can be obtained only in parameterized ways, faithfully 
taking into account ambiguities [521 [77]. At this dynamical stage, the construction of min- 
isuperspace operators currently proceeds by analogy with full operators, not by derivation. 

"'^^It is possible to construct flux operators from the volume operator, viewing the latter as some kind of 
basic operator However, for the present purposes one cannot substitute the volume for fluxes because 
no linear basic algebra would result for the definition of quantum representations and their averaging and 
reduction. 

^^Moreover, the volume operator usually used in homogeneous models, and also here, is a simplified 
version of the cubic SU(2)-invariant of the full theory. The assumed simplification of the much more 
complicated full spectrum does not follow from reduction but is put in by hand. When details of the 
eigenvalues are important, for instance when one uses the volume as the independent variable of difference 
equations, the simplified spectrum could lead to additional artefacts, not covered by the methods of this 
article. 
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3.6 Quantum-geometry corrections 

Using holonomies instead of curvature or connection components implies quantum-geometry 
corrections in the dynamics. There is a second type of effect, called inverse-triad correc- 
tion, which comes from the fact that an inverse of the densitized triad appears in the 
Hamiltonian constraint of gravity and in matter Hamiltonians, but flux operators have 
discrete spectra containing zero. No inverse flux operators exist, but the inverse densitized 
triad can be quantized to a densely defined operator using classical rewritings following 
[591 [78]. Ill Hamiltonian constraint, inverse-triad operators appear in the gravitational 
part (giving rise to differences of volume eigenvalues in ( 127]) and in matter Hamiltonians). 

Holonomy corrections are controlled by the parameters A/, or by the values 5/ chosen for 
a constraint operator. Inverse-triad operators entering ( l23l) via \p5ji,i/2{9K)[p5K,i/2{9K)~^ i 
are built using the same type of holonomies, and so their corrections refer to the same lat- 
tice scales 5/ as holonomy corrections. Both corrections are therefore linked to each other, 
and comparing the explicit forms of corrections allows one to estimate which one might be 
dominant in a given regime. 

Like holonomy corrections, the size of inverse-triad corrections depends on the values 
of 5i and requires a proper consideration of lattice structures. However, there is an addi- 
tional operator, the volume V , used crucially in the definition of inverse-triad operators 
as commutators. For this operator, the same question must be asked as for holonomies, 
namely what lattice scale it refers to. In a local lattice picture, as in the full theory, one 
should expect the relevant volume to be the one associated with a single lattice site or a 
spin-lattice vertex, just as the holonomies used correspond to single lattice links. However, 
incorporating the volume in this way is not as obvious as for holonomies, and so different 
versions have been considered, making use of macroscopic volumes [179] or even one asso- 
ciated with the artificial integration region V. In this subsection, we derive in detail the 
form of inverse-triad operators and the corrections they imply. To simplify commutator 
calculations involved, we will present the main equations for Abelian models and briefly 
comment on non-Abelian effects later. 

3.6.1 Local and non-local lattice operators 

Working with lattice spin-network states, one can define different flux operators which 
all give rise to the same flux when averaged to minisuperspace operators. This situation 
complicates constructions in pure minisuperspace models and has led to considerable con- 
fusion. Only relating models to the full theory, completing the kinematical reduction, can 
solve these issues. 

Local lattice operators: We begin with the local flux operator, able to show any inho- 
mogeneity realized in the lattice model: taken for a plaquette transversal to a 
surface and intersecting only one edge e^j starting at the vertex v. We choose 
the surface to be a square of coordinate area d-jd-K = ^j^kLjLk, so that we can 
view F^j as a quantization of the classical XjXkP^ (v), where p^{v) is an inhomoge- 
neous diagonal component making the homogeneous variables position dependent. 
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The conjugate variable A/Cj is quantized via local holonomies h^j = exp{i J cjds). 
(Recall our Abelian simplification in this subsection.) These local lattice operators 
satisfy the commutator algebra 

[K,i,Fy'j>] = -Sn-fil6^^y>6iyh^j . (28) 



Minisuperspace operators: If each surface used for local flux operators is centered at 
the intersection point with e^j, the union of all those that have the same /-coordinate 
as V form a surface stretching through the whole integration region V, without over- 
lap of non-zero measure. Including an average in the transversal direction, we can 
view the lattice sum = J\ff^ J2v according to ( 120|) as the flux quantizing the 
minisuperspace variable = LjLxP^ ■ Its conjugate variable c/ = LjCj in minisu- 
perspace is quantized by holonomies, hi = exp{ici) for an edge stretching through 
the whole integration region in direction X^. We have the commutator 

[hj,p^] = -87r-fel6fhi (29) 

for minisuperspace operators, correctly quantizing fll2p . 

Non-local operators: There is a version of operators between local lattice and minisu- 
perspace ones. We can average local lattice fluxes F^j over the lattice rather than 
sum as in p^ , or reduce the size of the minisuperspace flux p^ by dividing by the 
number of vertices in a surface, and define 

V 

This flux operator refers to the lattice spacing but, via averaging, includes all lattice 
sites in the integration region. We will call it the non-local flux operator. With a 
local holonomy, it obeys the commutator relation 



In (pT]) . the number A/1A/2A/3 of lattice sites in a region V replaces the coordinate 
volume Vq of ( !T2|) . At a technical level, l/NiN2N'^ comes about as the product of 
I/MjMk in the plane average (15U]) . and a factor of l/A/} because only one out of 
A/} lattice links along direction / provides a non-zero commutator [hyj,Fy/j] oc S^y 
according to f l28|) . 

It may seem questionable to use local holonomies and non-local fluxes within the same 
setting fl3T]) . but a consistent and closed algebra of basic operators is obtained in this way 
(provided the A/} are fixed). Whether such operators are meaningful physically is another 
question which we will soon discuss. For now, our motivation for looking at such a mix of 
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local and non-local operators is that it has been used (implicitly or explicitly) in several 
proposals to formulate inverse-triad corrections. 

Properties of basic operators in the different algebras can be translated into one another 

and are mutually consistent. If K/i/)) denotes non-local flux eigenstates with Fj| (///)) = 
SnjipfijKfij)) , the holonomy-fiux algebra (!3T!) determines the action h^jlui) = + 
I/A/1A/2A/3)) of local holonomies. Constant shifts of flux eigenvalues result for a fixed 
lattice. The quantized densitized-triad component = LjL^p^ is obtained from the 

averaged flux as = AfjAfxEi. Its eigenvalues change under the action of a basic holonomy 
operator by / = 8n-f i^XjAfKl^i ^ 87r74A/'jA/'x(/i/ + l/MiM2M^) = / + l/Mi. This 
dependence of the constant shift on Mi is consistent with the form exp{icj/Afi) of local 
holonomies exp(z£/C/), with cj = LjCj, ii/Lj = 1/A/}, and {ci,p'^} = 8n'yG6j . Notice the 

different behaviors of the non-local averaged flux Fj and the minisuperspace densitized- 
triad component p^ , corresponding by its factor of LjLk to the flux through a complete 
plane in the region V. 

The general form of the algebra of basic operators does not depend much on whether 

the local Fj, the non-local Fj or the minisuperspace p^ is used; the latter two differ from 
each other just by constant factors at the kinematical level (disregarding lattice refine- 
ment). However, Fj and p^ are much less local then F^j and therefore unsuitable for local 
expressions such as quantized Hamiltonians or inverse-triad corrections. 



3.6.2 Inverse-triad corrections 

The gravitational part of the Hamiltonian constraint contains a factor of e'^^^eabcEjE^/ ■>/ 1 det E 
in which one divides by the determinant of Ef, and similar terms occur in matter Hamil- 
tonians. Flux operators and the volume operator having zero in their discrete spectra, no 
densely defined inverse exists to quantize 1/ detE directly. Instead, one makes use of the 
classical identity [SU] 

27r7Gsgn(det E ^';^^' = {Al [ ./\d^\d'x} (32) 
^^\detE\ J 

and quantizes the Poisson bracket to a commutator of the form 

^ {h:'[K, V] - k[h-\ V]) = ^ {heVk' - k'vk) . (33) 

In a lattice model, holonomies refer to lattice links, or psj{gi) in their reduction. The 
volume operator is expressed via fluxes, and here local flux operators are used, given the 
local form of the classical Poisson bracket in (l32i) and of the commutator in (l33ll which 
depends only on vertex contributions to V lying on the edges used in hg. At this stage, 
minisuperspace models can easily become misleading because their most immediate flux 
operators p^ or p^ , proportional to Fj are non-local. The wrong form and size of inverse- 
triad effects then results. 
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We now present the detailed derivation of inverse-triad corrections based on local flux 
operators, as in [21 [80] , and then show how non-local versions differ. The simplified volume 
operator of Abelian models is = \Fv,iFv,2Fv,'i\^^'^ summed over all vertices of a 
spin-lattice state. In expressions such as ( l33l) . it suffices to look at contributions from all 
lattice-aligned hyj. A single such commutator is then 



Fvj\hv,i — h 



vJ 



F,,\h, 



v,I\"'v,I 



to be summed over all /. (Classically, the combination of holonomies and fluxes corresponds 
to \Fyj\~^/'^ ^J\F\^JF\k\. We have used an absolute value around the commutator as in 
( |23|) .) For Abelian holonomies it is easy to simplify the inverse-triad operator, making use 
of the commutator \hv,i-, F^j] = —Sii'yiphvj from ( l28l) and the reality condition h\ jhyj = 1. 

2 1 1/2 



Commuting holonomies past flux operators then gives hi j\Fyj\^^'^hyj 
and therefore 



l-v,! 



FvjFyxl 



(34) 



In strong quantum regimes, non-Abelian features should be relevant [81j and inverse- 
triad effect compete with holonomy and higher- curvature terms; however, the form ( IMl) 
still plays a characteristic role in effective actions The expression (IMl) is a good 
approximation in perturbative settings with F^ j > STT'jip, where it may be used to estimate 
qualitative effects or potential observational tests |82l [831 ED] • 

Since inverse-triad operators are local — commutators he[h~^, V] provide contributions 
only for vertices on e even if the volume operator for the full region V is used — their 
commutators refer to local F^j in = '^y\Fv,iFy^2Fv,3\^^'^, not to the minisuperspace 

operator or the non-local Fj. Inverse-triad corrections therefore depend on Fyj^STc-fip, 
where the Planckian addition can easily be a significant contribution to the eigenvalue or 
expectation value of F^j, the flux through an elementary lattice site0 Had we used the 

average Fj, the algebra would have led us to Fj ± Sn'yip / AfiAfjAfx , with corrections not 
only much suppressed by dividing by the large number of lattice sites but also depending 
on the size of the arbitrary region V choseni^i Such operators would be incorrect; they are 
based on the confusion of the correct average y/Fj with the non-local \/Fj. 



^^In fact, if Fyj is Planckian, with lattice spins near 1/2 for the fundamental representation, as often 
assumed, inverse-triad corrections are large. Geometry must be sufficiently excited above fundamental 
spins (some kind of ground state) for good semiclassical states to result. 

^^Sometimes, it is suggested to take a limit of Vb — )■ oo, or Afj — oo, viewing a finite Vb as a regulator. 
The procedure removes any V dependence and makes inverse-triad corrections disappear. However, as 
discussed in more detail in the next section, this reasoning is misguided: Vq is not a regulator because 
its value does not at all affect the classical theory. Classical models with different Vq produce the same 
physics, and so they should all be quantizable, without an effect of Vq. Moreover, the limit of A// — > oo is 
not consistent with the basic algebra of averaged operators. 
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3.6.3 Local quantum corrections 

We have distinguished three types of constructions for composite operators quantizing a 
symmetric model: the minisuperspace treatment using hj = exp{ici) and p"^ with algebra 
(!29|) . chimerical constructions with (local) link holonomies h^j but non-local fluxes Fj with 
algebra (I3T]) . and finally local lattice operators built from h^j and Fy^j with algebra ( 128|) . 

Local and minisuperspace treatments differ from each other by the order in which re- 
duction and composition of operators are done. In non-local models, as in traditional 
minisuperspace versions, one first postulates or derives the reduced basic operators hj and 

Fj (or p'^) and their algebra, and in a second step constructs composite operators of the 

form O non-ioca\{h I, Fj) from them by simple insertions, following analogous steps taken 
in the full theory. In local quantizations, one first constructs operators Oiocaiihyj, Ftv,j), 
adapting the full techniques to lattice states, and then restricts them to a quantized homo- 
geneous model. The second, local method is more complicated because it must deal with 
the reduction of non-basic, composite operators or their averaging. Tractable techniques 
exist only in rare cases, and therefore the main effects, for instance in the Hamiltonian 
constraint, are incorporated by parameterizations or mean-field techniques as in Sec. 13.5.31 
— an unsurprising feature given that local methods are analogous to a transition from 
microscopic Hamiltonians to tractable models of large-scale effects in condensed-matter 
physics. 

Despite technical difficulties, the local viewpoint has several clear advantages: it pro- 
duces the correct sizes of quantum corrections and naturally gives rise to lattice refinement. 
As already noted, the misrepresentation of quantum corrections in non-local models can 
easily be seen for inverse-triad operators, or the key ingredient O = Non-local 
operators make use of the averaged flux before taking the square root, quantizing O as 

(^non-local = |-^/|^'^^- A local quantization, by contrast, leads to (9iocai = iFvjl^^"^, the over- 
line now indicating restriction to homogeneous states after taking the square root. While 
linear combinations of the basic operators commute with averaging (Sec. I3.4.4p — produc- 
ing similar-looking basic algebras in local and non-local versions — non-linear combinations 
do not. In non-linear combinations of the basic operators, drastic deviations between local 
and non-local operators can therefore result, but only the local version correctly captures 
properties of the full theory in which no averaging is done. 

Local composite operators can be formulated only when the full lattice structure is 
taken into account, but after reduction they refer to reduced degrees of freedom as suitable 
for a quantization of a classically reduced symmetric model. Minisuperspace and non-local 
operators may be formally consistent without direct reference to a lattice, provided one 
chooses the length parameter ij = SiLj of holonomies in some way, for instance related 
to the Planck length or to full area eigenvalues, but at this stage the models become ad- 
hoc. Moreover, in spite of their formal consistency, non-local composite operators do not 
provide the correct form of corrections in operators that refer to fluxes or the volume, most 
importantly inverse-triad corrections. 

In addition to making inverse-triad corrections sizeable and interesting, the local treat- 
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ment taking into account inhomogeneity, has another imphcation regarding the physical 
evaluation of models. There is just one set of parameters, 5/ or equivalently A/} = 1/5/, 
which determines the magnitude of holonomy and inverse-triad corrections. It is not pos- 
sible to ignore one of the corrections and focus only on the other, unless one can show that 
a regime of interest leads to values of 5/ that make one correction dominate the other. In 
general, the two corrections are not strictly related to each other, because holonomy cor- 
rections are sensitive to the classical curvature scale relative to the Planck density, while 
inverse-triad corrections are sensitive to the local discreteness scale F^j relative to the 
Planck area, as seen in A detailed, state-dependent analysis, taking into account 

inhomogeneous quantum geometry, is required to estimate both corrections. 

4 Limitations of minisuperspace models 

Minisuperspace models of quantum cosmology never provide exact solutions to full quan- 
tum gravity. In some cases, deviations can be strong, for instance when unstable dynam- 
ics of neglected degrees of freedom (a classical property) enlarges the mismatch between 
symmetric and less-symmetric solutions, which at an initial time may have been the con- 
sequence only of a mild violation of uncertainty relations [84]. The discreteness of loop 
quantum cosmology shows a much larger class of minisuperspace limitations, for discrete- 
ness is not easily reconciled with homogeneity. As always, such limitations should be 
pointed out and discussed in detail, not to slander but to warn. 

At the level of states and basic operators, homogeneous wave functions can be derived 
in precise terms, using the distributional constructions of [19], as recalled in Sec. 13.4.21 
Intuitively, averaging a discrete state over a continuous symmetry group cannot result in 
a normalizable wave function in the original Hilbert space (or even a meaningful density 
matrix), but it is well-defined as a distribution. At this level, discreteness is not problematic 
and does not introduce ambiguities. At the dynamical stage, however, discrete space-time 
structures with possible refinement or (even if there is no refinement) reference to the local 
discreteness scale are more complicated and more ambiguous, as occasionally pointed out 
well before loop quantum cosmology was introduced [HSIISSIEZ]- Loop quantum cosmology 
has provided means to analyze such situations. 

4.1 Parameters 

In addition to phase-space variables, a strict minisuperspace model has only the parameter 
Vo = L1L2L3 to refer to, appearing in the symplectic structure ([1]). There is no analog of 
ii = XiLj = Lj/J\fi, the discrete lattice scales. And yet, the local quantum dynamics of 
the full theory, together with the quantum corrections it implies, depend on the parameters 
Xj via holonomies around loops used to quantize F*^, or along edges used in commutators 
with the volume operator to quantize inverse triads. 

It may not be obvious how exactly edge lengths XiLj enter quantum corrections, owing 
to a certain conceptual gap between the coordinate dependent Lj or and geometrical 



38 



aspects in this background-independent formulation, as well as quantization ambiguities. 
But quantum corrections certainly cannot depend on Vq, which is chosen at will (the 
coordinate size of a region used to reduce the symplectic structure) and knows nothing 
about the discrete scale. If Vq or any of its factors Lj appears in quantum corrections, an 
artificial dependence on coordinates and the chosen region results, as well as wrong sizes 
of quantum effects. 

In an ad-hoc way, [66] proposed to modify the strict minisuperspace dynamics in a 
fashion that successfully eliminates the Vo-dependence at least in holonomy corrections. 
(Inverse-triad corrections could not be represented meaningfully in this scheme.) As a 
consequence in isotropic models, the Hubble parameter rather than c = 70, appears in 
holonomies, and the discrete dynamics proceeds by constant steps of the volume V^a?, not 
of the densitized triad V^^^p. Heuristically, as in Sec. I3.5.5[ one can argue for this scheme 
by identifying geometrical areas a^£^, instead of coordinate areas with the Planck area 
when specifying the size of holonomy modifications. 

Holonomy corrections depending on the Hubble parameter a/ a rather than a have the 
advantage of being easily coordinate independent. If the Planck length — or a parameter 
close to it such as the smallest non-zero area eigenvalue of the full theory — is chosen as 
the discrete scale, modifications sin(£p7{)/£p oil-L result, independent of coordinates and of 
Vq. Holonomy corrections then refer simply to the curvature radius relative to the Planck 
length (or, via the Friedmann equation, the density scale relative to the Planck density), 
a parameter which can easily be estimated in regimes of interest. 

As a general scheme, however, the procedure suffers from several problems: 

1. If the volume is used as a basic variable, following [66], one introduces a sign choice 
by hand, allowing all real values for v = ±Voa^. (Otherwise, id/dV is not essentially 
self-adjoint, and exp(td/dV) not unitary. One would have to use the methods of 
affine quantum gravity [88] for acceptable quantizations, but it has not been shown 
that this can be compatible with the use of holonomies.) In fluxes, the sign appears 
automatically thanks to the orientation of triads [3^, and exp(td / dp) is unitary. 

2. The Planck length in sm{ip'H)/ip enters by an ad-hoc choice, which cannot be 
avoided because the quantization, still at the minisuperspace- level, does not have 
access to discrete structures. One may use the full area spectrum to guess what the 
scale might be, but such a procedure leaves open the question of what structure or 
eigenvalues a dynamical discrete state might give rise to. Moreover, the minisuper- 
space area or volume spectrum does not have a smallest non-zero eigenvalue. The 
spectrum, seen in ([8]), is discrete, with all eigenstates normalizable, but on the non- 
separable kinematical Hilbert space the spectrum still amounts to a continuous set 
of numbers as eigenvalues. One has to go slightly beyond minisuperspace models by 
referring to the full area spectrum, which does have a smallest non-zero eigenvalue, 
but still no reduction is performed. In this way, the scheme becomes improvised, 
heuristic, and ad-hoc. 

3. Going beyond strict minisuperspace quantizations is not a bad thing; in fact, it is 
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required for realistic modeling. However, improvised schemes do not provide jus- 
tifications for the detailed way in which one tries to go beyond. Moreover, while 
they give rise to meaningful results for holonomy corrections, inverse-triad correc- 
tions from (IMll are not modeled properly. These corrections depend on the ratio of 
the discreteness scale | (F) \ to the Planck area. If one assumes that the discreteness 
scale is exactly the Planck area, inverse-triad corrections would merely result in a 
constant factor, not affecting the dynamics much at first sight. But the factor dif- 
fers from one, and it has dynamical effects even if it changes just slightly. In other 
attempts of improvised quantizations, \{F)\ was related to macroscopic areas [79] . 
sometimes even involving Vq, for instance by using areas related to the size of the 
region V. These proposals ignore the fact that there is only one discrete structure 
that both holonomy and inverse-triad corrections can refer to, as well as their local 
nature. 

4. To counter inappropriate references to Vq in non-local quantizations, the parameter 
is sometimes treated gulator to be sent to infinity after quantization. Such a 
formal limit would undo all inverse-triad corrections, leaving only the classical inverse 
in dynamical equations. However, the limit does not exist at the level of operators 
— if it existed, it would result in an inverse of the triad operator, which is not 
densely defined. One can perform the limit at the level of the difference equation 
for wave functions, or in effective equations. But while this is formally possible, the 
overall quantization procedure would no longer be coherent. After all, Abelian loop 
quantum cosmology has difference equations for states because curvature is replaced 
by holonomies, resulting in a true modification sin{ii'H)/ii for "H. The limit — > 
or Xj = l/A/} —7- 0, which would send holonomy operators to derivatives by p, 
does not exist at the operator level. The Hamiltonian constraint is quantized to a 
difference operator with non-zero step-size. At the level of wave equations, acting 
with the operator on states in the triad representation, the limit does exist and 
produces a version of the Wheeler-DeWitt equation [89]. If one insists on removing 
the "regulator" Vq by sending it to infinity at the level of wave equations, one should 
also remove the true regulators ij in holonomies at the same level. The dynamics of 
loop quantum cosmology would then be no different from Wheeler-DeWitt dynamics. 

In the new homogeneous quantization of this article, the limit L/ — )■ cxd or Vq = 
L1L2L3 — )■ 00 is impossible at fixed 5/, because gi = exp{Li(l){Tj)) has no such limit. 
For edge lengths ij = 61L1 in psi,jj{gi) to remain finite, one would have to take the 
limit 5/ — > simultaneously with L/ — t- 00, but then the difference equation would 
become a differential equation, and loop quantum cosmology would, again, reduce 
to Wheeler-DeWitt quantum cosmology. Instead, one must be able to derive models 
for arbitrary values of Lj, such that observables are independent of them. 

5. The parameter Vq, in contrast to 61 in holonomy modifications, is not a regulator 
because it does not modify the classical theory. Classical models can be formulated 
with all finite choices of Vq, producing the same dynamics and observables. (Differ- 
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ent choices of Vq to some degree resemble different normalizations of the scale factor. 
Rescaling Vq is not a canonical transformation as it changes the symplectic structure 
(II]). Classically, this is not a problem, but one cannot expect a unitary transforma- 
tion at the quantum level, making the issue in quantum theory more complicated.) It 
should then be possible to formulate also quantum dynamics for all possible choices, 
or else Vq would acquire more physical meaning than it deserves. Another prob- 
lematic feature of regularization attempts is a possible topology dependence. If one 
looks at a model of closed spatial slices, for instance the FLRW model with positive 
curvature, one cannot send Vq to infinity. Instead, it may seem natural to use the 
full spatial coordinate volume as a distinguished value of Vq. (But again, one may 
equally well formulate the classical dynamics with different values of Vq, choosing 
different coordinates on the unit 3-sphere or integration regions smaller than the 
whole sphere.) It is sometimes argued that some effects, such as inverse-triad cor- 
rections, are meaningful or non-zero only with closed spatial topologies, but not for 
the flat, non-compactified FLRW model. Not surprisingly for a quantization based 
on non-local fluxes, quantum dynamics would then suffer from a strong violation of 
locality, depending on the global spatial topology even in its elementary changes. 
Such models in cosmology would also be hard to test empirically. One would have 
to know the spatial volume — and whether it is compact or not — before one can 
estimate quantum effects and make predictions. A more detailed discussion is given 
in[2]. 

The many conflicting comments that can be found in the literature of improvised quantiza- 
tions, for instance regarding the size of inverse-triad corrections, attest to the complicated 
and incomplete state of affairs in this scheme. Sometimes, a single paper may claim that 
inverse-triad corrections are too small to be significant, and at the same time can be 
changed at will by tuning the value of Vq. Although it is not always realized by all authors, 
such conflicting statements spell out limitations of pure minisuperspace models. 

4.2 Parameterizations 

The improvised scheme, by introducing "holonomies" as functions of the Hubble parameter 
rather than the connection component, mimics an a-dependent ij = XjLj oc 1/a = 1/ ^/\p\. 
As a or p changes and the universe expands or contracts, the lattice spacing evolves. 
Although there is no explicit creation of new vertices, the number Af of lattice sites must 
change, for a fixed Vq with changing A/ implies an evolving = I/A1A2A3. 

Deriving a precise functional form for \i{pi) would require one to formulate a corre- 
spondence between full discrete dynamics and reduced minisuperspace dynamics, including 
projections of evolved states onto the space of symmetric states. Lacking such complicated 
constructions, one can use phenomenological input to restrict possible forms of A/(p^), 
or \i{p'^) in diagonal anisotropic models, for instance different exponents of power- laws 
A(p) oc \p\^ in isotropic models with a real parameter x [35]. If A is constant (a; = 0, 
corresponding to ^39j), the discreteness scale would be magnified by cosmic expansion. 
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presumably making it noticeable in observations. Since no discreteness has been seen, 
X = or values close to it are ruled out. An improvised scheme amounts to x = —1/2, 
with constant discreteness scale, and is compatible with observations. However, a con- 
stant discreteness scale is not in agreement with full constraint operators changing vertex 
structures and local volume values, not just the number of vertices. On average over many 
individual actions of the Hamiltonian constraint and on large scales, cosmic minisuperspace 
dynamics may be close to x = —1/2 as in the toy model presented in Sec. 13.5.31 but this 
value cannot be realized precisely. 

The improvised scheme is compatible with most cases of cosmic evolution, but it has 
problems with black-hole models [90]. By its construction using geometrical areas of the 
region V, the scheme relates the number of vertices to the total volume of spatial regions. 
Near the horizon in homogeneous coordinates of the Schwarzschild interior, the spatial vol- 
ume shrinks, making the number of lattice sites small. However, the regime is supposed to 
be semiclassical for large black-hole mass, which is in conflict with a small number of lattice 
sites, implying noticeable discreteness. The analysis of different models — cosmological 
ones and those for black holes — shows that there cannot be a single universal power-law 
exponent for Xi{p'^) in all regimes. Discrete quantum dynamics and refinement behav- 
ior, just as the underlying state, depend on the regime analyzed. The role of coordinate 
choices hints at another important issue, namely how much the condition of covariance 
and anomaly freedom restricts possible refinement schemes. This question remains largely 
unexplored owing to the complicated nature of the quantum constraint algebra, but see 
[32] for an interesting cosmological example that suggests restrictions, also pointing at a 
value near x = —1/2. 

4.3 Reduction 

Lattice refinement in difference and effective equations refers to state parameters, most 
importantly Xj, depending on a geometrical variable such as the total volume V. One may 
view the appearance of V as an internal time, on which the evolving state depends. A 
possible procedure of implementing such a dependence, as alluded to in Sec. I3.5.3[ would 
be to write a full state as a superposition J2v '^v contributions ipv belonging to some 
fixed volume eigenvalue V. One would decompose a dynamical state as an expansion in 
eigenstates of the internal-time operator, such as the volume. Although the procedure 
would be difficult owing to the complicated volume spectrum and the fact that one would 
have to solve for a dynamical state first, it is in line with standard treatments of internal 
time. The states appearing in the decomposition of ipv in the spin-network basis then show 
what discrete structures are realized at a given volume V, and the spacing as well as the 
number of vertices might certainly change as one moves from one V to the next. 

After decomposing a dynamical state as '^yipv, one would still have to adapt it to 
near-homogeneous geometries, that is implement the projection back to symmetric states. 
Additional state-dependent parameters may arise, all to be modeled by suitable functions 
\i{Pj), the only parameters that survive with exact homogeneity. Such functions would 
then be inserted in dynamical equations of reduced models, for instance in difference equa- 
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tions of Abelian models. 

Notice that one must know some features of the full evolution of a state before defining 
the reduced Hamiltonian constraint, which in a second step can be used to evolve a reduced 
state. Since the Hamiltonian constraint is one of the operators to be averaged for reduction, 
evolution and reduction are not independent processes in the construction of models. As a 
consequence, reduced Hamiltonian constraints are state-dependent, even more so than the 
full constraint operator with its state-dependent regularization of j59] . 

If for a precise reduction we must know how to evolve a full state, why do we not 
work with the full evolved state rather than its reductions? The advantage of reduced 
models is that they offer additional approximation schemes, for instance in the derivation 
of observables or of effective equations. However, reduced models can never provide ex- 
act predictions — if their predictions were exact, one would not be dealing with a reduced 
model. It does not make much sense to derive physical quantities, for instance bounce den- 
sities, in exact terms within minisuperspace models because the models themselves are not 
exact. Only general effects, such as quantum hyperbolicity, the presence of bounces under 
certain conditions, or qualitative low-curvature corrections may be meaningful predictions, 
but not specific values of some parameters related to the discreteness scale. 

4.4 Spin-foam cosmology? 

Spin-foam cosmology [9ll [921193] attempts to enlist spin-foam techniques to address quantum- 
cosmological questions by embedding a simple structure with finitely many edges in a spa- 
tial (or space-time) manifold E. Such a map is clearly different from the (mini-)superspace 
embedding — )■ 5 used for classical reductions, or a map a: Tihom — ^ Anhom of state 
spaces used for quantum reductions. Looking back at our discussion in Sec. Ej the question 
therefore arises what kind of construction spin-foam cosmology can provide. 

As usually emphasized in this context, spin- foam cosmology aims at a description of 
quantum cosmological space-times without making use of reductions, rather describing 
physics in a full theory of quantum gravity in which inhomogeneous modes are still present 
and quantum fluctuate. According to the classification in Sec. Ej such a non-reduction 
scheme could only be selection or projection, but there is certainly no control over the full 
non-symmetric solution space, let alone the averaging problem, within spin-foam cosmol- 
ogy. Spin-foam cosmology does not fall within our classification of reduction schemes. 

In fact, it is not clear in which sense — or if at all — spin- foam cosmology describes 
symmetric models. It is true that inhomogeneous modes have not been truncated but 
remain present and may fluctuate, potentially a feature that would allow one to go beyond 
reduced models (see also [9ll|95] in the canonical setting). However, spin- foam cosmology 
at present lacks conditions that would ensure inhomogeneous modes to be sufficiently small 
for the models to be considered symmetric, not just in their fluctuations but even in their 
expectation values. The graphs used in spin-foam cosmology refer to finitely many degrees 
of freedom, often related heuristically to the number of degrees of freedom of homogeneous 
minisuperspaces. However, counting degrees of freedom is not enough to ensure that a 
model is good. One might simply define a finite-dimensional "minisuperspace" by picking 
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some point xq in space-time and considering only the metric components gab{xo) as de- 
grees of freedom, a model which would be meaningless because of its dependence on the 
space-time gauge. By using spatial embeddings, of graphs instead of points xq, spin-foam 
cosmology is in danger of producing models close to the one just sketched. (Indeed, the 
status of covariance remains unclear in spin-foam cosmology as well as full spin foams.) 

What is missing in this context is a well-defined analog of the map a for states used 
in loop quantum cosmology. Only such an object could tell whether the correct degrees 
of freedom have been captured. Another question, related to the topics of this article, is 
how spin foams Abelianize. The final equations often produced in this context resemble 
difference equations of Abelian loop quantum cosmology, even though the starting point 
has SU(2) degrees of freedom. No clear Abalianization step has been provided. Finally, 
working with fixed graphs embedded in space, spin- foam cosmology has not given rise to 
refinement models. 

5 Evaluation of models of quantum cosmology 

So far, we have discussed the construction of reduced and other models. Their limited 
nature regarding the dynamics requires care also, and especially, when they are evaluated 
for physical predictions. In addition, there are caveats which apply to any construction in 
quantum gravity, and so to model systems as well. To guarantee that models, obtained 
in a sufficiently parameterized way to ensure their genericness, can indeed be evaluated 
reliably, one must use evaluation methods or approximations that do not bring in hidden 
assumptions about the general form of effects. We end this article with a brief exposition 
of the main features; see also the review [96] . 

5.1 Problem of time 

As with all totally constrained systems, the problem of time [971 [981 199] , plays an important 
role to be addressed when one tries to evaluate theories. In models of quantum cosmology, 
deparameterization is an oft-used method to obtain an evolution picture with corresponding 
observables [lOOj . possibly related to physical predictions. However, deparameterizability 
is not generic and requires specific, unrealistic matter choices. If deparameterization is 
used to extract predictions, one would have to make sure that results do not depend 
on assumptions made about matter fields or internal times 1^ Deparameterization often 
ignores potentially important off-shell effects and, in the presence of inhomogeneity, has no 
direct means of testing whether a quantum version of the contracted Bianchi identity holds 
— an important condition for the consistency of Einstein's equation or any modification or 
quantization thereof. Moreover, the scheme presupposes that some notion of time remains 

"'^^Models in which physical Hilbert spaces are derived based on deparameterization are often caUed 
"complete quantizations." However, as a quantization of a space-time theory, such constructions can be 
considered complete only when one has shown that results do not depend on the choice of internal time. 
No such demonstration has been given in the models proposed so far. See also the discussion in [lOlj . 
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valid even in strong quantum regimes, and is therefore blind to a possible signature change, 
as it occurs in loop quantum cosmology at high density [9]. 

Although deparameterization provides established methods to derive physical Hilbert 
spaces, the independence of results on deparameterization choices, a form of space-time 
anomaly problem, has never been demonstrated so far. However, there are effective proce- 
dures in which different such choices are realized as gauge transformations at least semiclas- 
sically |102[ I103[ 1104] . At the semiclassical level, which is sufficient for most questions of 
direct physical interest with some observational contact, the problem of time can therefore 
be overcome. 

5.2 Effective equations 

The classical /imit presented in Sec. 4.3 of |40], based on [89], has often been misinterpreted 
as a semiclassical approximation. Accordingly, important quantum corrections have been 
missed in many analyses: While some /i-terms are kept in £p-related holonomy corrections, 
quantum back-reaction terms of the same order are dropped. 

Quantum back-reaction is captured by effective actions or effective equations of motion. 
Canonical effective equations [731 1105] available at the semiclassical level are derived from 
expectation values {C) of constraint operators in a general class of semiclassical states 
[Til 11061 1107] . Properties of these states are derived in this scheme along with quantum- 
corrected solutions for expectation values and need not be presupposed, a feature that helps 
ensure genericness. Expectation values of the Hamiltonian constraint are subject to three 
types of quantum corrections: holonomy and inverse-triad corrections from the specific form 
of quantum geometry realized in loop quantum gravity, and, as in all interacting quantum 
theories, quantum back-reaction of fluctuations and higher moments on expectation values. 

Quantum back-reaction is a dynamical feature and more complicated to derive than 
holonomy and inverse-triad corrections. But it is an important effect, giving rise to higher 
time derivatives in effective equations, one ingredient of higher-curvature terms in effective 
actions of quantum gravity. Just as holonomy corrections, quantum back-reaction effects 
are therefore expected to depend on the curvature scale — both types of corrections should 
be of similar sizes. (However, holonomy corrections are not identical to or the sole contrib- 
utor to curvature corrections.) In models in which quantum back- reaction terms have not 
been computed, holonomy corrections cannot be seen as reliable: higher-curvature terms 
may undo any holonomy effect. In particular, if the whole series of holonomy modifications, 
such as sin(£/C/)/£/ instead of c/, is used, an inconsistent approximation is made. One must 
view iCj) / 1 1 = ci{l + 0{ijcj)) as a series expansion, in which each term competes with 
one of similar size but of higher-derivative form, for instance cj compared with c/. (When 
inhomogeneity is included, also higher spatial derivatives result from holonomies of inte- 
grated y4Jj, expressed as a derivative expansion to collect correction terms of similar sizes, 
but they have rarely been included in calculations so far.) If quantum back-reaction terms 
are unknown, strong holonomy effects, for instance around the maximum of the sine as 
used in bounce models, are not reliable. 

The only models in which quantum back-reaction terms are known at present are har- 
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monic ones, entirely free of quantum back-reaction [108[ [72] . In all other cases, the high- 
curvature regime remains uncontrolled. (If one evolves whole wave functions, for instance 
numerically, quantum back-reaction is implicitly included. Several non-harmonic models 
have been analyzed in this way ^8 11109] . However, this requires specific assumptions about 
an initial state to be evolved, usually assumed to be closely related to a Gaussian, a sin- 
gle or at most two-parameter family within the infinite-dimensional (semiclassical) state 
space. Since it is difficult to guess what a quantum state at high curvature should look 
like, current numerical methods are not generic enough for robust predictions.) 

5.3 High- curvature regime 

Ignorance of the high-curvature regime would seem to preclude conclusions about singular- 
ities. Fortunately, there are generic statements that can be made, even though they do not 
give detailed intuitive pictures such as bounces. The most general statement in Abelian 
models (mini- [HI |39l [Ml |52] or midisuperspaces [1101 1111] ) is quantum hyperbohcity |7], 
stating that the wave function can be extended uniquely beyond geometries that would 
classically be singular. (In non-Abelian models, where difference equations are lacking, 
quantum hyperbohcity is realized by the annihilation of singular states by the Hamilto- 
nian constraint.) The statement applies to all solutions of the difference equation, and is 
therefore independent of the structure of physical Hilbert spaces or wave-function normal- 
ization. It is free of limitations related to deparameterization. The conclusions obtained 
with quantum hyperbohcity apply not just in homogeneous but also in midisuperspace 
models. Moreover, they extend to generic inhomogeneity using a new scenario |112] rem- 
iniscent of that due to BKL |113j . Thanks to inverse-triad corrections, spatial derivatives 
in effective actions are suppressed close to spatial singularities, in a more generic way than 
classical BKL statements would suggest. 

The same methods show that large holonomy corrections imply signature change [9]: 
at strong curvature, there is only a quantum version of 4-dimensional Euclidean space, no 
space-time if holonomy corrections are dominant. This effective picture may seem to be 
inconsistent with an "internal-time evolution equation" for states such as ( 1271) . which takes 
a wave function through the small-volume regime in which holonomy effects are strong. 
However, in these high-density regimes the right-hand side of ([271) is very large. Solutions 
of the difference equation then no longer oscillate because the oscillation length would 
become much smaller than allowed by the discrete spacing. Such regimes are classically 
forbidden, as can be seen heuristically by the fact that the bounds provided by holonomy 
functions sm.{liCi) j would be surpassed at high density. The classically forbidden regime 
of the difference equation indicates that space-time and evolution as we know them are no 
longer realized, just as an effective action with signature change would imply. Tunneling 
from negative to positive or nj (or between two opposite orientations) and a Euclidean 
phase are two sides of the same mechanism, seen only when inhomogeneity is properly 
taken into account. The signature-change scenario is in fact in agreement with a tunneling 
regime in quantum hyperbohcity, also providing interesting parallels with the pictures of 
|114[ 1115] . as already anticipated in p.l6j . 
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6 Conclusions 



Much work remains to be done to establish reduced models of loop quantum gravity as 
well-defined and controlled approximations, and as reliable sources of detailed predictions 
in high- curvature regimes. As discussed in Sec. [21 an analogous problem must be faced even 
classically in relating non-symmetric geometries to symmetric ones: the averaging problem. 
A complete understanding of quantum minisuperspace models as approximations of the 
full theory, even if a reduction mechanism is included, can be obtained only when the 
classical averaging problem is better understood. Lacking a general solution, no approach 
to quantum gravity is yet able to produce a complete derivation of reduced models. 

Nevertheless, with sufficient care one can make progress and render it at least likely 
that all crucial effects of the full theory are captured. That non-Abelian effects should play 
some role in loop quantum cosmology and require caution has been emphasized quite some 
time ago regarding specific properties of inverse-triad corrections [HI] as well as general 
properties of homogeneous models [6l], but it has not often been realized. The present 
article offers several new observations and constructions to this end: We have pointed 
out that most considerations made so far in loop quantum cosmology suffer from Abelian 
artefacts, related to the use of function spaces on the Bohr compactification of the real line, 
following 110]. To correct this oversight, a new quantization of homogeneous connections 
is developed in Sec. [31 which starts from non-Abelian models and takes into account the 
complete structure of invariant connections. The resulting Hilbert-space representation, 
when restricted to Abelian variables, is related to Bohr-quantized models by an isometric 
*-algebra morphism, but one that is not unitary or bijective. We lose information when we 
map states to the Bohr Hilbert space, corresponding to the edge-spin degeneracy inherent 
in previous models. 

The edge-spin degeneracy of holonomies is removed by the new quantization in non- 
Abelian and Abelian models, giving a better handle on lattice structures and the relation 
to the full theory. We have strictly related basic operators in the full theory and in 
models, and showed how quantum effects in composite operators can be captured by local 
quantizations. The averaging required can lead to unexpected features — as seen in detail 
for fiux operators — which one would not endeavor to implement in a pure minisuperspace 
quantization without being prompted by the relationship with the full theory. Several 
implications have been demonstrated, especially regarding lattice refinement and the form 
and sizes of quantum-geometry corrections, most importantly those due to inverse triads. 

We emphasize that our constructions started with the realization of deficiencies in 
current quantizations based solely on Abelian models; seeing lattice refinement or local 
features was not the main aim but nevertheless resulted as an unavoidable consequence. 
Sees. [H and [5] have provided cautionary remarks, detailing the current incomplete status 
of the field and providing some guidelines for evaluations and the approach to physicality. 
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